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ABSTRACT 


When  a  beam  of  extraordinarily  polarized  light  is  directed  through  a  poled,  single 
crystal  of  BaTiOj  it  fans  in  the  direction  of  the  optic  axis.  This  thesis  research 
investigates  the  origin  of  this  fanning,  or  asymmetric  self-defocusing,  in  barium 
titanate.  Experimental  observations  and  theoretical  simulation  suggest  that  this 
phenomenon  is  due  to  near  forward  stimulated  photorefractive  scattering. 
Additionally,  observations  of  unique  variegated  beam  patterns  during  phase 
conjugation  are  reported,  including  internally  pumped  rings  and  motion  invariant 
patterns. 
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I.  INTRODUCTION 


Dunng  the  last  decade  the  photorefractive  crystal  barium  titanate  (BaTiO,)  has 

- - - - ^ 

received  a  large  amount  of  attention  within  the  scientific  commumty.  Barium  titanate 

has  long  been  of  interest  to  acousticians  and  solid  state  physicists  because  of  it's 

photoacoustic  and  piezoelectric  properties,  but  only  recently  have  optical  physicists 

become  very  interested  in  this  crystal.  The  cause  of  this  sudden  and  intense  interest 

was  the  discovery  of  photorefraction  and  optical  phase  conjugation.  Originally  found 

in  liquids  and  gases  (Ref.  1],  optical  phase  conjugation  has  the  ability  to  reconstruct  a 

distorted  wavefront  (sometimes  loosely  referred  to  as  time  reversal).  Barium  titanate  is 

one  of  many  photorefractive  crystals  which  can  exhibit  phase  conjugation. 

There  exist  a  number  of  intriguing  optical  properties  of  single  crystals  of  barium 
titanate,  but  by  far  the  most  interesting  is  its  ability  to  phase-conjugate  at  continuous 
wave  low  powers,  with  no  external  apparatus.  This  effect,  known  as  self-pumped  phase 
conjugation  (SPPC),  is  the  subject  of  extensive  current  investigation.  The  goal  of  this 
research  is  to  determine  the  origin  of  self  pumped  phase  conjugation  in  BaTiO^. 

Currently  there  exist  two  explanations  of  the  cause  of  SPPC  in  BaTiOj.  One 
theory,  developed  by  Feinberg  and  his  students  at  UCLA  (Ref  2),  attributes  the 
production  of  the  phase  conjugate  reflection  to  the  photorefractive  effect  via 
Degenerative  Four  Wave  Mixing  (DFWM).  According  to  this  theory,  an  asymmetric 
change  in  the  index  of  refraction  across  the  path  of  an  incident  beam  causes  the  beam 
to  bend.  This  bending  eventually  leads  to  a  geometry  which,  through  two  internal 
reflections  (often  termed  a  CAT  comer),  allows  the  beam  to  intersect  itself,  leading  to 
DFWM.  The  second  theory,  originating  with  Lam  and  his  co-workers  at  Hughes 
Research  Laboratory  (Ref  3],  attributes  SPPC  to  Simulated  Photorefractive  Scattering 
(SPS),  or  Two  Wave  Mixing,  similar  to  that  exhibited  in  gases  and  liquids  with  two 
frequency  disparate  input  beams  of  very  high  powers.  With  the  goal  of  proving  one  or 
the  other  of  these  theories,  this  thesis  begins  with  a  general  background  and  theory  of 
phase  conjugation  as  it  applies  to  BaTiOj  (Chapter  II).  Chapter  III  describes  PSIM 
(Photorefraaive  SIMulation),  a  simulation  of  the  phenomenon  of  beam  fanning  due  to 
a  change  in  index  of  refraction  within  the  incident  beam  path  arising  from  the 
photorefractive  effect.  This  simulation  enables  the  prediction  of  the  change  in  index  of 


refraction  and  the  index  profile  necessary  to  account  for  SPPC  with  only  the 
photorefractive  effect  taken  into  account. 

Chapter  IV  introduces  experiments  to  determine  the  effect  of  incident  beam  size 
on  various  crystals  of  BaTiOj.  Chapter  V  summarizes  the  results  and  analyzes  the 
data  presented.  The  conclusion  is  reached  that  Stimulated  Photorefractive  Scatten.ng  ;s 
responsible  for  the  process  of  beam  fanning  which,  leads  to  self  pumped  phase 
conjugation  in  BaTiO^.  The  actual  origin  of  the  phase  conjugate  beam  is  attributed  to 
both  SPS  and  DFWM. 

Additionally,  Appendix  A  contains  a  summary  of  the  work  I  did  at  Los  Alamos 
National  Laboratory  under  the  direction  of  Dr.  R.  A.  Fisher  in  February  and  March, 
1986.  This  work  describes  and  catalogs  internal  beam  production  in  barium  titanate 
during  SPPC.  Some  internal  beam  patterns  not  previously  reported  are  introduced  and 
explanations  for  these  phenomena  are  proposed.  The  results  of  the  experimental  work 
in  the  body  of  this  thesis  have  been  submitted  for  publication  to  Optics  Letters.  An 
expanded  version  of  Appendix  A  was  presented  at  the  XIV  International  Quantum 
Electronics  Conference  (1986),  and  has  been  submitted  for  publication  in  the  Journal  of 
the  Optical  Society  of  America  B  in  collaboration  with  R.  A  Fisher  and  A.  V.  Nowak. 


II.  BACKGROUND  AND  THEORY 


A.  INTRODUCTION 

Optical  phase  conjugation,  first  discovered  in  1972  [Ref.  1],  relies  on  a  specific 
nonlinear  propeny  of  a  substance  to  create  the  formation  of  the  complex  conjugate  of 
an  incident  wave.  The  reflection  involved  is  not  the  familiar  specular  reflection  seen  in 
ever>'day  life,  but  a  reflection  that  retraces  the  path  of  the  incident  beam  exactly, 
regardless  of  the  direction  or  angle  of  incidence.  The  exact  retracing  of  the  incident 
beam  path  and  the  formation  of  the  complex  conjugate  of  the  incident  wavefront  are 
the  two  unique  and  most  important  effects  that  make  phase  conjugation  both 
interesting  and  useful.  In  barium  titanate  this  phase  conjugate  beam  is  directly  a  result 
of  the  photorefractive  effect. 

B.  ORIGINS  OF  THE  PHASE-CONJUGATE  BEAM 

The  origin  of  the  phase  conjugate  beam  in  photorefractors  is  light  which  is  Bragg 
scattered  by  an  index  grating  formed  by  the  light  in  the  material.  This  process  is 
known  as  wave  mixing.  There  exists  both  two-wave  mixing  (scattering)  and  four-wave 
mixing  in  barium  titanate.  Both  of  these  processes  will  be  treated  theoretically  in 
general  terms;  the  specifics  of  wave  mixing  in  BaTiOj  will  be  treated  later. 

To  understand  the  origins  of  phase  conjugation  one  must  begin  with  the  wave 
equation  for  a  homogeneous,  isotropic,  time-independent,  nonmagnetic  and 
nonconducting  dielectric.  The  polarization  term  is  divided  into  its  linear  (Pj)  and 
nonlinear  (Pj^j)  parts,  and  the  linear  portion  is  absorbed  into  the  linear  permittivity  of 
the  medium.  £.  This  leaves  the  wave  equation  as; 

V^E  -  f2.1) 

cf  dr 

For  ease  of  manipulation  an  appro.ximation  to  this  equation  is  made  that  reduces 
the  second  order  nonlinear  wave  equation  to  a  first  order  equation.*  This 
approximation  is  derived  by  Fisher  and  Yariv  [Ref  5:  pp.  9-11]  and  for  it  to  be  valid, 
the  envelope  of  the  pulse  must  not  change  appreciably  during  an  optical  period,  hence 

*Some  calculations  can  be  made  without  this  approximation  but  they  are  quite 
tedious  (e.g.  [Ref  4]). 


the  term  Slowly  Varying  Envelope  Approximation.  In  making  the  approximation  it  is 
necessary  to  assume  that  the  electromagnetic  waves  under  consideration  are  plane 
waves,  and  also  to  extract  only  the  portion  of  the  nonlinear  polarization  that 
propagates  with  the  same  frequency  and  wave  vector  as  the  modulation  function  of  the 
electric  field.  This,  referred  to  as  the  phase  matched  portion  of  the  polarization,  is 
given  by 

Pni  ®  exp ± i((ot -  kz),  (2.2) 

in  the  wave  given  by 


E  s  E(z,t)  exp ±  i(<ot  —  kz).  (2.3) 

Here  the  boldface  letters  (E  and  P)  are  the  fields  and  the  conventional  letters  (E  and  P) 
are  the  field  envelopes,  which  vary  little  in  an  optical  period  and  in  an  optical 
wavelength.  All  other  terms  are  ignored  because  they  are  not  phase  matched  and 
therefore  cannot  couple  to  the  electric  field.  The  Slowly  Varying  Envelope 
Approximation,  or  SVEA,  is  then  invoked.  That  is: 
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Fisher  and  Yariv  then  show  that  the  wave  equation  with  SVEA  becomes  a  first-order 
equation  relating  the  envelope  functions  E  and  P: 
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This  resuit  is  the  SVEA  '.vave  equation  for  a  plane  wave  traveling  in  the  -  z  direction, 
which  shows  that,  given  a  nonlinear  polarization,  an  electric  field  can  be  produced. 

This  development  can  be  generalized  to  the  case  of  many  waves  present  within  a  single 
medium,  resulting  in  the  origin  of  the  nonlinear  processes  that  can  lead  to  phase 
conjugation. 

The  polarization  term  can  often  be  written  as  a  power  series  expansion: 
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P  -  EX(E)  -  X(‘^E  -t-  X^2)e2  ^  x<^)E^  -t-  •  , 


(2.6) 


where  X{E)  is  the  susceptibility  of  the  material  in  which  the  wave  is  traveling.  The 
expansion  of  the  susceptibility  in  a  power  series  separates  the  effects  of  the  nonlinear 
matenal  into  distinct  groups  that  categorize  nonlinear  elTects  [Ref.  5:  p.  Ijj: 

•  X^*^  -  These  are  the  linear  properties  that  are  the  subject  of  classical  optics  and 
are  not  directly  responsible  for  nonlinear  effects. 

•  X^^^  -  The  second  order  effects  are  in  general  called  three-wave  mixing  and 

include:  second-harmonic  generation,  optical  rectification,  parametnc  mixing 
and  the  Pockels  effect. 

•  X^^^  -  Third  order  effeas  include  some  of  the  most  popular  effects  responsible 
for  phase  conjugation  such  as;  third-harmonic  generation,  nondegenerate  four- 
wave  mixing,  Raman  scattering,  ac  Kerr  effect  (degenerate  four-wave  mixing), 
dc  Kerr  effect,  Brillouin  scattering,  two  photon  absorption,  dc-induced 
harmonic  generation  and  Stimulated  Photorefractive  Scattering. 

The  two  methods  of  creatmg  a  phase-conjugate  reflection  that  are  directly  applicable  to 

the  experimentation  described  here  arc  scattering  and  degenerate  four  wave  mixing. 

1.  Four- Wave  Mixing 

In  1977  Hellwarth  (Ref  6)  showed  that  it  was  possible  to  generate  a  time- 
reversed  replica  of  any  monochromatic-beam  wave  pattern.  The  method  employed  the 
interaction  of  an  incident  beam  with  counter-propagating  pump  waves  m  a 
homogeneous,  transparent,  nonlinear  medium  (see  Figure  2.1).  This  method,  called 
four-wave  mixing,  has  become  one  of  the  most  popular  methods  for  phase  conjugation 
because  of  the  almost  universal  availability  of  the  necessary  equipment  and  proper 
media. 

The  description  of  wave  mixing  begms  with  the  nonlinear  polarization  derived 
from  the  energy  function.  As  shown  by  Pepper  and  Yariv  [Ref  5:  pp.  26-331.  isolating 
only  the  third  order  term  gives 
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.  vCi.s.ucr  i'iCti  regenerate  iav..r-.va'. e  nu.xnig  vDl  'A  vi  .  xnere 
ojj  -  Wj  =  (Op  =  (0^  represents  frequencies  of  the  four  beams  The  subscripts  I  and  2 
refer  to  two  strong,  precisely  counterpropagating  pump  beams,  the  subscript  p  refers  to 
the  weak  probe  beam,  and  the  subscript  c  refers  to  the  phase-conjugate  output  \va\e. 
It  IS  assumed  that  the  following  conditions  are  met: 


Figure  2.1  Schematic  of  the  process  of  four-wave  iruxing  [Ref.  5:  p.4211. 

a)  irttensc  (nondcpletcd)  counterpropagating  planewave  pump  beams  of  equal 
intensity  (kj,  k,), 

b)  a  weak  probe  beam  (kp), 

c)  both  the  intense  pump  beams  and  the  weak  probe  beam  are  at  the  same 
frequency  to. 

Since  the  two  pump  waves  are  very  strong  compared  to  the  probe  beam  they  form  a 
standing  wave.  To  obtain  the  steady-state  solutions,  we  can  set 


(2.8) 


So.  assuiraiis:  :ae  SVE.\  and  considering  only  the  phase  matched  portion  of  the  phase 

conjugate  wave,  then  from  Equation  2.5 


Where 


(2.16) 


E^.(L)  is  the  amplitude  of  the  conjugate  wave  at  its  origin  (the  phase  conjugate  mirror), 
and  Ep(0)  is  the  amplitude  of  the  probe  wave  at  its  origin. 

In  practice,  the  amplitude  of  the  phase  conjugate  wave  at  its  origin  is  zero  (i.e. 
E^(L)  =  0),  so  that  at  the  input  point  (z=0)  of  the  probe  beam 

E,(0)  =  i  — tan(|k|l)Ep*(0).  (2.17) 


Equation  2.17  shows  that  at  z=0  (the  origin  of  the  probe  beam)  the  reflected  field  is 
proportional  to  the  complex  conjugate  of  the  incident  field.  Therefore  there  exists  a 
phase  conjugate  reflection. 

2.  Scattering 

Phase  conjugation  by  stimulated  scattering  was  first  reported  by  Zel'dovich, 
Popovichev,  Ragul'ski  and  FaizuUov  in  1972  [Ref.  1].  As  stated  by  Hellwarth 
(Ref.  5:  p.  171]; 

In  essence,  the  nonlinear  polarization  density,  mediated  by  some  driven  (Raman-, 
Brillouin-,  etc.)  active  excitation  of  the  optical  medium,  couples  normal-mode 
solutions  of  the  linear  Maxwell  equation  so  as  to  create  a  new  set  of  normal 
modes  in  backscattering  with,  in  effect,  complex  propagation  constants 
representing  index  change  and  gain. 

For  example  Stimulated  Brillouin  Scattering,  or  SBS,  is  a  direct  result  of  the 
electro strictive  effect  seen  in  some  materials.  This  effect  is  seen  when  the  polarizability 
of  a  material  is  a  function  of  pressure.  In  these  materials  it  is  possible  for  light  to  be 
scattered  by  a  sound  wave  allowing  the  coupling  of  a  pair  of  light  waves  to  a  pressure 
wave.  Another  type  of  scattenng  sometimes  responsible  for  phase  conjugation  is 
Stimulated  Raman  Scattering  or  SRS.  SRS  is  the  scattering  of  light  from  molecular 
vibrations  or  rotations  in  a  gas  or  liquid,  or  longitudinal-optical  phonons  in  a  solid.  In 
SRS  light  waves  are  coupled  because  of  polarizability  changes  with  a  molecular 
coordinate.  Of  particular  importance  to  this  research  is  Stimulated  Photorefractive 


Scattering  in  barium  titanate.  In  SPS  it  is  the  photorefractive  effect  which  is 
responsible  for  the  production  of  a  periodic  index  of  refraction  grating;  this  grating 
then  Bragg  scatters  light  in  the  backward  direction  which  comprises  the  phase 
conjugate  beam. 

The  mathematics  involved  in  all  three  types  of  scattering  is  similar  and  can  be 
treated  together  .  To  this  end  I  outline  the  work  of  Hellwarth  (Ref.  5:  pp.  177-180]  and 
discuss  the  existence  of  the  phase  conjugate  wave  assuming  scalar  or  symmetric 
scattering  with  a  monochromatic,  multimode,  incident  wave  in  a  wave  guide.  The 
results  are  largely  independent  of  the  nature  of  the  wave  guide  provided  that  a  specific 
TEM  mode  that  propagates  in  the  forward  direction  also  propagates  in  the  backward 
direction. 

A  wave  E^(r)  may  be  expressed  as: 

Ei(r)  *  f  exp[i(k^z  -  pz)/2] ,  (2. 18) 

where  p  is  the  attenuation  coefficient,  6^  is  the  normalized  transverse-mode  pattern, 
and  the  propagation  constants  may  vary  with  wave  angular  frequency,  that  is  k^^  ■ 

Uv)- 

The  electric  field  of  the  backscattered  wave  may  then  be  written  as 
[Ref.  5;  p.l80] 

^2  ”  exp[-ik„z  +  pz/2  -  yz/2] ,  (2.19) 


where  k^^  *  kjj(<o),  and  y  is  used  to  balance  any  nonlinear  terms  from  the  fields 
interaction  with  the  medium.  If  y*0  in  the  above  equation,  then  this  is  the  solution 
to  the  linear  Maxwell  equations  for  a  wave  in  a  wave  guide. 

If  in  Equation  2.19  y*0.  then  there  is  one  configuration  that  allows  the  y  of 
one  of  these  modes  to  have  a  positive  real  part  significantly  greater  than  that  of  any 
other.  This  means  that  a  single  wave  (reflection)  may  dominate  the  backscattered 
waves  if  the  necessary  conditions  are  met.  The  desired  result  is  that 


Ejfr)  »  nEj  (r), 


(2.20) 
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where  n  is  a  constant.  Substituting  Equations  2.18  and  2.19  into  Equation  2.20  gives 
the  polarization  density  necessary  to  balance  the  extra  terms  in  Maxwell's  equations 
that  come  from  the  last  term  in  Equation  2.19  .  According  to  Hellwarth  [Ref.  9]  the 
necessary  conditions  are  satisfied  if 


(2.21) 


where 


mijn 


||^(1/L)  exp[AkLz]dz, 


(2.22) 


and 

Ak  -  kj(0)  +  ki(v)-kj(v).k„«0).  (2.23) 

In  Equation  2.21  I  have  assumed  that  the  scattering  considered  is  the  pure  scalar  or 
symmetric  type,  reducing  the  fourth  rank  susceptibility  tensor  to  a  scalar  quantity. 
This  quantity,  and  all  other  constants,  I  have  absorbed  into  n. 

The  necessary  conditions  may  be  determined  for  special  cases,  either  by  exact 
solution  or  penurbation  methods  developed  by  Hellwarth.  These  solutions  indicate 
that  a  backscattered  solution  that  is  nearly  phase  conjugate  predominates  because  it 
has  gain  nearly  twice  that  of  any  other  mode,  provided  that  [Ref  5:  p.  177]: 

a.  the  interaction  length  is  not  too  long, 

b.  the  Stokes  shift  (v  -  <o)  is  not  too  large, 

c.  the  number  of  guide  modes  falls  within  a  certain  large  range,  and 

d.  the  total  number  of  guided  modes  is  not  too  large. 

The  solutions  that  led  to  these  conditions,  presented  by  Hellwarth  in  197S  [Ref.  91, 
untbrtunateiy  Jo  nor  take  into  .iccount  competing  nonlinear  eifects,  and  therefore  ire 
good  only  to  the  extent  that  they  predict  the  major  effects  observed  in  experimentation. 
It  is  observed  that  the  backscattered  wave  is  not  always  an  exact  phase  conjugate  of 
the  input  wave,  possibly  due  to  some  of  these  competing  effects.  These  analyses  apply 
equally  as  well  to  SPS,  SBS  and  SRS  [Ref  3). 


C.  DISTORTION  CORRECTION  BY  PHASE  CONJUGATION 

Although  not  immediately  applicable  to  this  investigation,  the  ability  of  a  phase 
conjugator  to  correct  for  phase  distortions  is  one  of  the  primary  aspects  driving 
research  in  this  area.  As  a  practical  point,  this  ability  is  an  easy  way  to  determine  that 
any  given  observed  beam  is  a  phase  conjugate  of  another. 

It  was  shown  above  that,  at  least  for  the  two  cases  considered,  a  wave  can  be 
generated  propagating  directionally  opposite  to  an  input  wave,  with  it's  envelope 
function  proportional  to  the  complex  conjugate  of  the  input  beam's  envelope  function. 
This  phase  conjugate  wave  may  be  considered  a  separate  wave  that  must  satisfy  the 
same  wave  equation  as  the  input  wave.  This  means  that  if  the  wavefront  of  the  input 
wave  is  distorted  in  some  way,  then  the  phase  conjugate  wave  will  propagate 
backwards  along  the  same  path  and  its  envelope  function  will  remain  everywhere  the 
complex  conjugate  of  the  input  wave's;  that  is,  its  wave  fronts  uHll  coincide  with  the 
input  wave  fronts  at  every  point.  Thus  a  distorted  wave  front  will  be  restored  to  its 

original  shape  as  the  phase  conjugate  wave  traverses  the  distorting  medium  in  the 

• 

reverse  direction.  The  proof  is  quite  simple,  as  shown  by  Fisher  and  Yariv  [Ref.  5:  p. 
17],  and  is  accomplished  by  complex  conjugating  the  wave  equation  for  the  forward 
going  wave.  The  result  is  the  equation  for  the  conjugate  wave  propagating  in  the 
opposite  direction.  So  a  wave  traveling  in  the  opposite  direction  of  the  incident  wave, 
with  an  envelope  function  the  complex  conjugate  of  such  a  wave,  satisfies  the  same 
wave  equation  as  the  incident  wave.  This  means  that  the  backward  going  wave  will 
remain  everywhere  the  complex  conjugate  of  the  incident  wave,  and  therefore  a 
distorted  plane  wave  will  be  restored  to  its  original  configuration  when  it  is  passed  back 
along  the  same  distorting  path.  This  distortion  correction  capability  is  the  most  visible 
and  fascinating  effect  attributable  to  phase  conjugation  at  the  present  time. 

An  incomplete  list  of  some  applications  of  phase  conjugation  with  emphasis  on 
those  currently  using  barium  titanate  as  the  phase-conjugator  is  [Ref.  5,10,11]: 

a.  Holography 

b.  Distortion  correction 

c.  Laser  osciiiators 

d.  Interferometry 

e.  Laser  beamsteering 

f  Associative  memory  investigations 

g.  Image  amplification 

h.  Laser  coupling 


i.  Image  addition,  subtraction  and  enhancement 

j.  Optical  bistability 

k.  Automatic  pointing  and  tracking. 

D.  PHOTOREFRACTORS  AND  FOUR  WAVE  MIXING 

In  photorefractive  materials  the  index  of  refraction  changes  with  incident  optical 
energy.  This  effect  was  first  discovered  in  the  1960's  [Ref.  12,13]  and  has  since  led  to 
the  development  of  optical  recording  sensitivities  comparable  to  that  of  the  silver  halide 
emulsions  used  in  photography  [Ref.  14].  The  process  involved  in  the  effect  is 
sununarized  by  Feinberg  [Ref.  5:  p.  418]  as  follows: 

a.  Light  causes  charge  to  migrate  and  separate  in  a  crystalline  material. 

b.  The  separation  of  charge  produces  a  strong  electrostatic  field. 

c.  The  electrostatic  field  causes  a  change  in  the  refractive  index  of  the  crystal  by 
the  linear  electro-optic  effect.^ 

In  the  photorefractive  effect  the  important  faaor  in  the  change  in  index  of 
refraction  is  not  the  intensity  of  the  incident  light,  but  the  total  energy  of.the  incident 
radiation;  the  power  only  determines  the  speed  of  the  reaction.  Also  because  of  the 
type  of  effect  concerned,  the  actual  result  is  dependent  on  the  relative  intensity  of  the 
incident  radiation. 

Consider  a  thin  beam  of  light  incident  on  a  photorefractive  crystal  that  is 
otherwise  dark.  The  light  will  liberate  charge  carriers  from  the  lighted  area  of  the 
crystal  and  these  carriers  will  then  settle  back  into  a  place  within  the  crystal  structure. 
If  the  carriers  settle  back  in  the  area  of  the  crystal  that  is  lighted  they  will  again  be 
liberated.  If,  however,  the  carriers  land  in  a  dark  area  of  the  crystal,  they  will  stay 
where  they  land  due  to  the  small  dark  conductivity  of  the  crystal.  Eventually  there  will 
be  few  charge  carriers  in  the  lighted  area  and  many  in  the  dark  area.  This 
disproportionate  distribution  creates  an  electric  field  between  the  dark  area  and  the 
lighted  area  that  creates  a  change  in  the  index  of  refraction  of  the  crystal  through  the 
linear  electro-optic  effect. 

Tlus  process  of  migration  has  been  aptly  described  bv  R.  A.  Fisher  as  'he 
flyswatter  eiTect  [Ref.  16].  It  is  as  u  the  charge  carriers  were  ilies  at  a  picnic  table  and 
the  light  was  a  flyswatter  suspended  from  above,  slowly  swinging  across  the  table  top. 

^The  linear  electro-optic  effect  is  the  change  in  the  index  of  refraction  caused  by 
(and  proportional  to)  an  applied  dc  electric  field.  This  effect  is  seen  in  crystals  that  do 
not  exhibit  inversion  symmetry,  such  as  barium  titanate  in  the  tetragonal  phase 
[Ref  15:  p.  275]. 


Each  time  a  fly  lands  on  the  table  within  the  area  of  the  flyswatter's  reach,  he  will 
immediately  fly  up  as  the  swatter  returns  toward  him.  If,  however,  the  fly  lands  on  the 
table  in  an  area  away  from  the  swatter  he  will  be  left  alone  and  will  stay  on  the  potato 
salad.  In  time  therefore,  there  will  be  no  flies  on  the  table  in  the  area  where  the 
flyswatter  can  reach  and  many  where  it  cannot.  This  process  effectively  pumps  the 
flies  to  regions  on  the  table  which  are  away  from  flyswatter. 

The  model  of  charge  transport  described  above  is  known  as  the  hopping  model 
and  was  first  proposed  by  Chen  to  explain  the  change  in  refractive  indices  in  LiNbOj 
and  LiTaOj  [Ref.  17].  There  also  exists  a  diffusion  model  proposed  by  Amodei 
[Ref.  18]  that  may  dominate  when  the  spacing  between  light  and  dark  areas  within  the 
crystal  is  less  than  one  micrometer. 

From  the  flyswatter  analogy  it  is  clear  that  what  is  important  is  the  optical 
intensity  relative  to  the  ambient.  This  entire  explanation  assumes  the  existence  of 
charge  carriers  in  a  crystalline  material.  In  most  materials  exhibiting  the 
photorefractive  effect  the  origin  of  the  charges  is  unknown,  but  it  is  assumed  that  they 
inhabit  low-lying  traps  formed  by  impurity  sites  in  the  crystal.  For  example,  in 
BaTiOj,  Fe^'*’  and  Fe^"*"  impurities  may  create  the  charge  carriers  as  shown  by  Klein 
and  Schwartz  [Ref  19],  It  has  also  been  proposed,  and  there  is  extensive  evidence  to 
support  the  theory,  that  oxygen  vacancies  are  the  donors  and  barium  vacancies  are  the 
acceptors  [Ref  20].  Whatever  the  dominant  species  is,  it  is  assumed  that  there  are  a 
large  number  of  empty  receptor  sites  in  the  material  available  to  the  liberated  charge 
carriers. 

If  instead  of  a  single  beam,  two  beams  of  identical  wavelength  and  equal 
intensities  (one  called  a  pump  wave  and  one  called  a  probe  wave)  intersect  within  a 
photorefractive  crystal,  the  interference  pattern  they  form  will  cause  a  periodic  optical 
field  to  be  set  up  in  the  region  of  intersection,  resulting  in  a  periodic  charge  separation. 
Through  the  electro-optic  effect,  a  periodic  variation  in  index  of  refraction  will  result. 
This  index  grating  is  then  able  to  Bragg  scatter  a  pump  wave  in  the  opposite  direction 
to  produce  a  phase  conjugate  beam.  The  formation  of  the  index  grating  is  shown  in 
7'gure  l.Z,  adapted  from  Feinberg  [Ref  5:  p.  422). 

It  is  important  to  note  that  photorefractive  phenomena  are  not  merely  four- wave 
mixing  phenomena.  Four-wave  mixing  is  a  consequence  of  an  index  change  at  some 
given  point  at  some  given  time  (x,t).  Photorefractive  phenomena  however,  are 
dependent  on  the  history  of  the  medium  and  are  not  spatially  local  effects.  That  is,  the 
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Figure  2.2  Index  grating  formation  due  to  the  photorefractive  effect  (Ref.  5). 

effects  at  (x,t)  are  a  consequence  of  actions  at  (.\-6x,  t-5t).  This  is  clearly  seen  in 
Figure  2.2. 

The  top  curve  of  Figure  2.2  shows  the  intensity  as  a  function  of  the  position 
along  the  path  of  the  beam.  The  second  curve  is  the  charge  density  due  to  this 
intensity  variation.  The  electric  field  intensity  is  presented  in  the  third  curve  (note  the 
important  V4  period  shift  relative  to  the  charge).  Finally  the  variation  of  the  actual 
index  of  refraction  within  the  cr>'stal  is  shown  in  the  last  curve. 

The  spatial  shift  between  the  change  in  index  of  refraction  and  the  variation  in 
intensity  due  to  the  interference  pattern  (one  fourth  of  the  grating  period)  is  an 
important  effect.  This  happens  because  the  peak  of  the  charge  distribution  is  a  center 
or'  rigiit-left  symmetry,  and  thereibrc  mere  can  be  ao  electric  Held  at  that  point.  This 
shift  leads  to  the  ability  to  transfer  power  from  one  beam  to  another,  a  process  known 
as  two-beam  coupling.  Two-beam  coupling,  first  reported  by  Staebler  and  Amodei 
[Ref  21],  is  due  to  the  scattering  of  one  beam  into  another  by  Bragg  diffraction  as  they 
pass  through  the  grating  they  have  mutually  produced.  The  result  is  that  one  beam 
experiences  gain  at  the  expense  of  the  other. 


One  interesting  effect  of  four-wave  mixing  in  photorefractors  is  seen  in  the  phase 
conjugate  mirror  with  continuous  wave  gain.  This  was  first  demonstrated  by  Feinberg 
in  1980  [Ref  22]  using  the  configuration  shown  in  Figure  2.3  . 
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Figure  2.3  Formation  of  an  optical  cavity  with  a  phase  conjugate  mirror  [Ref.  5|. 

In  this  configuration  a  phase  conjugate  mirror  (PCM)  and  an  ordinary  mirror  (Ml) 
can  form  an  optical  cavity  in  which  a  beam  will  oscillate.  One  of  the  more  dramatic 
demonstrations  of  this  type  of  oscillation  is  reported  by  Feinberg  (Ref  5:  p.  436]  in  a 
configuration  using  an  ordinary  kitchen  spatula  as  one  end  of  the  cavity.  Just  as  in 
any  conventional  laser  oscillator,  the  oscillations  can  begin  with  noise  (as  little  as  a 
single  photon)  that  makes  one  successful  round  trip  witliin  the  cavity.  With  each 
additional  round  trip  the  system  sees  gain  due  to  two-beam  coupling  and  in  a  sliort 
period  of  time  a  visible  CW  beam  may  be  observed  between  the  mirror  and  the  phasc- 
conjugator. 

E.  STIMULATED  PHOTOREFRACTIVE  SCATTERING 

The  second  theory  which  accounts  for  phase  conjugation  in  BaTiO,  was 
jrigmaily  proposed  'oy  Lam  :n  198.5  [Ref.  jj.  Lam  proposed  diat  tne  origm  )i'  ;;ie 
phase  conjugate  beam  is  two  wave  mixing  or  scattering.  This  process  ol'  SPS  begins 
with  random  inhomogeneities  in  the  photorefractive  medium  which  give  rise  to  noise 
photons  from  the  Rayleigh  scattering  of  the  input  beam.  In  the  words  of  Lam  [Ref  J]; 

The  noise  photons  propagating  along  with  [sic]  axis  of  the  input  beam  undergo 
parametric  scattering,  resulting  in  an  unidirectional  gain.  The  phase  conjugate 


wave  arises  from  the  component  of  the  scattered  wave  which  experiences  the 
maximum  gain. 

Later  the  same  year  Chang  and  Hellwarth  conclusively  demonstrated  SPS  in  BaTiOj 
(Ref.  23].  To  do  this  they  had  to  index  match  the  surfaces  of  the  crystal  to  make 
accessable  angles  which  would  not  produce  the  characteristic  comer  reflection  (i.e.  the 
CAT  corner)  producing  DFWM. 

Later,  Valley  [Ref.  24]  defined  SPS  as  the  selective  amplification  through  the 
photorefractive  effect  of  optical  radiation  scattered  by  medium  inhomogeneities.  Thus, 
SPS  is  a  two-wave  mixing  process  in  which  the  incident  beam  interferes  with  scattered 
light  to  create  the  necessary  grating  through  the  photorefractive  effect.  As  soon  as  the 
gain  for  an  individual  noise  beam  exceeds  its  loss,  then  that  beam  will  become 
dominant,  increasing  the  grating  strength  and  hence  the  backscattered  beam  gain.  This 
process  will  continue  until  a  steady  state  value  is  reached. 

The  mathematical  equations,  as  derived  by  Lam,  are  of  the  identical  form  to 
those  encountered  in  stimulated  scattering  in  a  waveguide  (presented  above),  so  the 
analysis  of  phase  conjugation  by  SPS  has  already  been  addressed. 

F.  BARIUM  TITANATE 

1.  The  Crystal 

The  experiments  reported  here  were  performed  with  one  of  the  more  popular 
photorefractive  crystals  capable  of  creating  a  phase  conjugate  reflection.  Above  the 
Curie  temperature  of  120®C,  BaTiOj  is  a  cubic  crystal  and  therefore  centre  symmetrical 
and  non-photorefractive.  Between  5®C  and  120®C  however,  the  crystal  is  in  the  polar 
phase  and  has  a  tetragonal  symmetry.  The  axis  of  fourfold  rotation  is  the  optical  axis 
or  c-axis.  In  the  tetragonal  phase  barium  titanate  is  ferroelectric,  photorefractive, 
photoconducting,  birefringent,  electro- optic,  pyroelectric,  piezoelectric,  photoacoustic 
and  (slightly)  electro strictive.  The  many  possible  effects  make  the  crystal  itself  very 
difficult  to  understand,  since  often  the  results  of  the  various  effects  cannot  be 
separated.  All  results  reported  here  involve  BaTiO^  in  the  tetragonal  phase. 

Since  3a  FiO,  ;s  chermcaily  and  mec.hanicaily  stable  at  room  temperature,  it  is 
used  extensively  in  many  capacities  and  is  therefore  well  studied.  Jona  and  Shirane 
[Ref  25;  pp.  108-215]  report  extensively  on  the  structure,  properties  and  the  effects  of 
temperature,  stress  and  electric  fields  on  barium  titanate  in  the  tetragonal  phase. 
Although  at  the  time  of  printing  the  photorefractive  capabilities  were  unknown,  their 
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data  represent  a  comprehensive  picture  of  the  crystal  and  forms  the  basis  for  analysis 
of  the  solid  state  aspects  of  barium  titanate.  More  recent  but  less  comprehensive  data, 
emphasizing  the  optical  aspects  of  BaTiOj,  have  been  compiled  by  Wemple,  et.  al. 
[Ref  26]. 

When  used  as  a  photorefractor,  it  is  necessary  that  a  single  crystal  of  BaTiO^ 
be  used  and  that  the  crystal  be  poled  into  a  single  domain.  Poling  is  accomplished  by 
cooling  the  crystal  through  the  Curie  point  while  an  electric  field  on  the  order  of  10^ 
V/cm  is  applied  across  the  optical  axis.  This  poling  of  the  crystal  results  in  an  electro¬ 
optic  tensor  which  has  a  very  large  [42]  element.  Since  the  photorefractive  effect  uses 
the  linear  electro-optic  effect,  it  is  important  to  maximize  this  [42]  element.  This 
maximization  is  accomplished  by  using  extraordinarily  polarized  light  incident  at  an 
angle  between  zero  and  90°  to  the  c-axis. 

2.  Beam  Fanning 

When  a  beam  of  extraordinarily  polarized  light  is  incident  on  a  poled,  single 
crystal  of  barium  titanate,  it  is  seen  to  asymmetrically  defocus  in  the  direction  of  the  c- 
axis  prior  to  the  onset  of  phase  conjugation.  This  process,  called  beam  fanning,  was 
first  discovered  by  Feinberg  in  1981  [Ref  27].  Beam  fanning  is  clearly  visible  in  Figure 
2.4  (the  optic-axis  is  designated  by  the  arrow).  This  fanning  can  be  explained  by  either 
an  asymmetric  change  in  index  of  refraction  within  the  incident  beam  path,  or  by 
stimulated  scattering. 

a.  Deflection  by  Photorefraction 

Beam  fanning  was  proposed  to  be  due  to  an  asymmetric  change  in  the 
index  of  refraction  caused  by  the  photorefractive  effect.  This  process  is  represented  in 
the  curves  of  Figure  2.5  (from  Feinberg  [Ref  27]  ).  The  top  curve  shows  the  intensity 
curve  of  a  Gaussian  beam.  The  second  curve  shows  the  net  charge  due  to 
displacement  of  charge  carriers  by  the  incident  light.  The  third  curve  shows  the 
induced  electric  field,  and  the  last  curve  shows  the  asymmetric  change  in  the  index  of 
refraction.  It  is  assumed  that  at  the  center,  and  most  intense  portion  of  the  beam,  the 
gradient  of '.he  index  of  'efr.action  is  greatest.  Therefore  *he  center  portion  of  the  beam 
creates  tne  cbserv-ed  nfect.  and  tne  fanning  ;s  inuuced  beginning  at  that  point.  In 
other  words,  beam  fanning  begins  in  the  center  of  the  beam.  This  is  one  of  the  major 
discriminating  factors  and  will  be  referred  to  later.^  The  second  major  point  to  consider 


^Note  that  even  if  the  intensity  profile  of  the  incident  beam  is  not  smooth,  the 
top  curve  of  Figure  2.5  represents  the  envelope  that  the  beam  intensity  must  fall 
within. 


Figure  2.4  BaTiOj  exhibiting  SPPC. 
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about  this  theory  is  that  as  the  diameter  of  the  beam  decreases  the  electric  field 
increases,  since  the  linear  distance  between  charge  distributions  decreases.  Since  the 
photorefractive  effect  is  proportional  to  the  electric  field,  then  beam  fanning  must 
increase  with  a  decrease  in  beam  diameter.  This  is  clearly  seen  by  examining  the 
equations  describing  the  electric  field  and  the  subsequent  refractive  index  change. 

The  electrostatic  field  induced  by  an  incident  beam  in  a  photorefractor  is 
given  by  [Ref  27]; 


-2/criT(.\x  +  v-y)  -  ,  , 

E(.\)  =  — 2 - -  exp[(x^  +  y^),'(o  2], 


(2.24) 


•.vhere  ■;  :s  '.he  charge  of  'he  cr.arge  carrier,  -'vg  is  Soitzmann  s  constant,  T  :s  '.no 
temperature  of  the  crystal,  and  <Oq  is  the  incident  beam  diameter.  The  change  in  index 
of  refraction  in  a  crystal  of  BaTiO^  for  an  extraordinary  ray  due  to  this  field  is 
approximately  given  by 
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Figure  2.5  The  process  of  beam  fanning  due  to  photorefraction 
for  a  single  incident  beam  (Ref  27) 

A  COS0  ,  , 

-(E(x)«yl  — — — (2n  -n  -r^,sin-0). 

2n(0)  ^ 


(2.2.5) 


where  n(0)  *  (n^*  cos0  +  n^^  sin0)'*  “  ,  0  is  the  angle  the  incident  beam  makes  with 
the  non  c-a\is  of  the  crystal,  n^  is  the  index  of  refraction  for  an  ordinary  ray.  n^  is  the 
index  of  refraction  of  an  extraordinary  ray,  and  all  terms  that  do  not  include  the  r^j 
element  of  the  electro-optic  tensor  have  been  dropped,  since  its  value  exceeds  all  other 
elements  by  an  order  of  magnitude. 

Inspection  of  these  equations  clearly  shows  that  the  change  in  index  of 
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tl.  tl-  SPS 


In  19S5  Lam  proposed  that  beam  fanning  is  due  to  an  elTcctivc  nonlinear 
index  of  refraction  which  causes  the  crsstal  to  behave  like  a  diverging  lens  (Ref  5) 
Ibis  nonlinear  index  of  refraction  i  ^  is  due  to  an  intensity  dependent  fiequency 
shift  of  the  incident  beam  due  to  the  bulk  photovoltaic  elfect.  Recently  \'al!e\ 


(Ref.  24)  has  suggested  that  beam  fanning  is  actually  near  forw'ard  .scattering  due  to 
SPS,  in  an  extension  of  Lam's  theory’. 

As  in  other  types  of  scattering,  SPS  is  the  selective  amplification  of 
radiation  which  is  scattered  by  inhomogeneities  in  the  medium.  The  amplification  is 
accomplished  through  the  photorefractive  elfect.  This  is  schematically  diagrammed  in 
Figure  2.6. 
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Figure  2.6  The  process  of  SPS  [Ref  24]. 

In  a  crystal  where  one  charge  carrier  dominates,  the  gain-length  product 
per  absorbed  energy  per  volume  at  turn  on  is  given  by  [Ref  24] 


/:refp«?firE^lgnr 

2nJcc^(l  +  E^/E^)|  ’ 


(2.26) 
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k  =  27t  X,  e  is  the  carrier  charge,  T  is  the  space  charge  field  and  refractive  index  grating 
decay  rate,  c  is  the  dielectric  constant  (which  is  anisotropic),  and 
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where  /cg  is  Boltzmann's  constant,  ^  is  the  carrier  mobility  (which  is  also  anisotropic) 
and  T  is  the  temperature  of  the  crystal.  Finally,  lg|j-  is  the  effective  interaction  length 
which  Valley  models  with 


Igfj-  “  L{1  -  exp[-<o^/(L  sine)]}, 


(2.28) 


where  L  is  the  length  of  the  crystal  in  the  direction  bisecting  the  scattered  and  incident 


beam,  co^  is  the  incident  beam  diameter,  and 
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(2.29) 


The  important  item  to  note  is  that  the  gain-length  product  is  proportional 
to  the  effective  interaction  length  given  by  Equation  2.28  .  Therefore  the  amount  of 
scattering  due  to  SPS  (both  forward  and  backward)  is  exponentially  dependent  on  the 
incident  beam  diameter  in  such  a  way  that  an  increase  in  beam  diameter  will  increase 
the  scattering,  and  a  decrease  in  beam  diameter  will  decrease  the  scattering;  this  is 
opposite  from  that  expected  by  the  theory  discussed  in  Section  F.2.a  above.  Also, 
since  SPS  is  dependent  on  the  interaction  of  random  noise  photons  with  the  incident 
beam,  one  would  expect  that  beam  fanning  (near  forward-scattering)  could  commence 
at  the  periphery  of  the  beam  in  addition  to  the  center.  This  scattering  of  energy  should 
be  more  obvious  at  the  edges  of  the  incident  beam  since  the  intensity  is  lower  there. 
Also,  any  light  scattered  toward  the  c-axis  from  inside  the  beam  would  take  power 
from  the  periphery  through  two  beam  coupling,  thus  adding  to  the  loss  of  intensity 
toward  the  beam  edge.  Considering  these  points,  if  the  intensity  profile  of  the  incident 
beam  were  to  be  monitored  during  the  farming  process,  one  should  expect  to  see  the 
intensity  diminish  from  the  periphery  toward  the  center.  Again  this  is  opposite  the 
effect  expected  from  the  prior  theory. 

3.  Self  Pumped  Phase  Conjugation 

Seam  fanning  is  intricately  involved  with  self-pumped  phase  conjugation 
(SPPC).  SPPC  in  BaTiOj  was  first  discovered  by  Feinberg  soon  after  the  discovery'  of 
beam  fanning  [Ref.  2).  According  to  Feinberg,  the  process  of  beam  fanning  (due  to  an 
asymmetric  change  in  index  of  refraction  within  the  beam  path)  may  result  in  a  portion 
of  the  incident  beam  being  bent  in  such  a  way  that  it  is  internally  reflected  off  the 
comer  of  the  crystal  (see  Figure  2.4).  This  reflected  beam  then  intersects  the  incident 


beam  causing  a  steady  state  interference  pattern.  The  light  and  dark  areas  of  the 
interference  pattern  cause  an  index  grating  (through  the  photorefractive  effect)  which 
then  Bragg  scatters  the  incident  beam  in  the  process  of  DFWM,  creating  a  phase 
conjugate  beam.  Lam  alternatively  theorized  that  the  origin  of  the  phase  conjugate 
beam  is  two-wave  mixing  (SPS)  as  described  above,  and  that  auxiliary  beams  were  not 
essential  to  the  process. 

Of  interest  also  is  the  time  for  phase  conjugation  to  begin  in  barium  titanate. 
The  time  to  produce  self-pumped  phase  conjugation  can  vary  from  a  few  seconds  to 
tens  of  minutes  depending  on  the  intensity  of  the  incident  beam  and  the  configuration 
used.  Because  of  the  typically  large  time  for  the  onset  of  phase  conjugation  in  BaTiOj, 
there  is  extensive  work  presently  underway  to  try  to  decrease  it.  Some  work  on 
changing  crystals  after  growth  has  been  reported  [Ref.  20],  but  the  present  emphasis  is 
on  the  doping  of  the  crystals  during  growth  [Ref.  29]. 

The  exact  processes  involved  in  beam  fanning  and  phase  conjugation  are  still 
unresolved  in  the  literature.  In  the  following  chapters  I  will  examine  these  theories  as 
they  apply  to  observed  effects  and,  since  in  the  case  of  beam  fanning  there  are 
mutually  exclusive  indications,  I  will  show  that  beam  fanning  must  be  due  solely  to 
near  forward  SPS. 


III.  PSIM:  A  SIMULATION  OF  THE  PHOTOREFRACTIVE  EFFECT 

IN  BARIUM  TITANATE 


A.  INTRODUCTION 

To  investigate  the  photorefractive  effect  in  barium  titanate,  I  designed  a 
computer  simulation  which  plotted  the  internal  beam  pattern  within  the  crystal.  This 
simulation,  entitled  PSIM  (Photorefractive  SIMulation),  only  accounts  for  the 
nonlinear  photorefractive  effect,  not  for  SPS  or  other  competing  nonlinear  effects.  In 
so  doing,  I  isolated  the  photorefractive  effect  from  other  effects  and  determined  the 
necessary  change  in  index  of  refraction  across  the  beam  path  to  account  for  beam 
fanning  which  is  typically  observed  in  the  laboratory.  The  source  code  for  PSIM  is 
found  in  Appendix  B. 

B.  THEORY 

As  discussed  above,  the  photorefractive  effect  may  be  solely  responsible  for  the 
observance  of  beam  fanning,  which  is  in  turn  responsible  for  self-pumped  phase 
conjugation  [Ref  2].  The  change  in  index  of  refraction  of  a  given  crystal.  An,  may  be 
given  by  the  equation  (Ref  14] 

An(z)= -(V2)n\{|'E(z),  (3.1) 

where  n  is  the  index  of  refraction,  r^f^^  is  the  effeaive  Pockels  coefficient  and  E(z)  is  the 
space  charge  electric  field. 

Values  for  r^jj- range  generally  in  the  picometers  per  volt  range.  Of  seven  crystals 
used  by  Klein  and  Schwartz  [Ref  19],  r^jp  ranged  from  4.2  x  10'^^  to  12.0  x  lO'^^ 
meters  per  volt.  Although  these  values  are  very  small,  the  space  charge  electric  field 
may  be  very  high,  on  the  order  of  10^  volts  per  meter  [Ref  5:  p.  418].  Using  these 
values  it  becomes  obvious  that  the  change  in  index  of  refraction  due  to  the 
photorefractive  elfect  in  barium  titanate  may  be  as  high  as  10"*  .  Though  not  an  exact 
number  for  any  given  crystal,  this  may  be  accepted  as  close  to  an  upper  bound  for  An 
(denoted  here  as  An^^^j^j^)  for  most  of  the  crystals  used.  Gunter  [Ref  11:  p.  226]  gives 
An^j^ax  ^  by  the  possible  space  charge  field. 


From  the  above  simple  calculations  it  may  be  assumed  that  for  the 
photorcfractive  effea  to  be  solely  responsible  for  self-pumped  phase  conjugation,  the 
dimensions  of  any  given  crj'stal  used  must  be  such  that  a  change  in  index  of  refraction 
on  the  order  of  10'^  is  sufficient  to  cause  the  observed  beam  fanning.  However,  as 
previously  mentioned,  the  gradient  of  the  electric  field  and  hence  the  index  gradient  is 
the  important  quantity  relating  to  the  photorefractive  effect.  Therefore,  the  beam  size 
relative  to  the  length  of  the  crystal  along  the  direction  of  propagation  and  the  intensity 
profile  must  also  play  an  important  role  in  SPPC.  These  factors,  in  addition  to  the 
incident  angle,  were  eventually  incorporated  into  the  simulation  described  here.  The 
resulting  simulation  was  used  to  determine  the  minimum  change  in  index  of  refraction 
necessary  to  account  for  observed  beam  fanning  in  any  given  configuration. 

C.  THE  SIMULATION 
I.  Program  Logic 

The  simulation  is  written  in  HP  basic  and  designed  to  run  on  a  Hewlett 
Packard  300  model  computer  with  at  least  two  megabytes  of  random  access  memory. 
This  two  dimensional  model  takes  into  account  only  the  change  in  index  of  refraction 
across  the  incident  beam  path  due  to  the  photorefractive  effea.  The  output  is  a 
graphic  representation  of  the  eventual  steady  state  beam  geometry. 

The  model  uses  a  brute  force  approach  to  calculate  the  observed  effects  which, 
although  time  consuming  and  memory  intensive,  uses  a  simple  algorithm.  The 
program  uses  only  Snell's  law  and  Huygen's  principle  to  calculate  the  internal  beam 
pattern.  Input  into  the  program  is  the  range  of  the  change  of  index  of  refraction 
across  the  input  beam.  An,  the  crystal  and  beam  size  (in  arbitrary  units),  the  point  of 
entrance  of  the  beam  into  the  crystal,  and  the  angle  of  incidence  of  the  incident  beam 
relative  to  the  normal.  The  program  computes  the  mean  index  of  refraaion  between 
two  input  extremes  which  then  corresponds  to  the  normal  index  of  refraction  (n  2.4). 
The  important  factor  in  the  program  is  the  difference  between  the  two  extremes, 
AOj^ax  •  Logically  it  is  this  entity  which  drives  the  photorefractive  effect.  Note  that 
the  simulation  is  not  concerned  with  how  :s  created. 

Once  the  necessary  parameters  are  entered,  the  internal  angle  of  the  mcident 
beam  is  calculated  and  each  point  on  the  wavefront  is  then  assigned  an  index  of 
refraction  value  that  it  can  impart  to  the  crystal  at  its  location.  These  points  are  then 
grouped  into  pieces  of  the  wavefront  (one  to  five  separate  pieces  are  normally  used). 


The  assigned  index  value  is  dependent  on  An^^^  and  an  assigned  profile  which  may  be 
varied  (more  will  be  said  about  this  profile  later).  The  indices  of  each  point  of  the 
wave  front  are  then  propagated  in  a  straight  line  through  the  crystal,  at  which  time  the 
value  of  the  index  of  each  point  in  the  wavefront  is  imparted  to  the  corresponding 
point  within  the  crystal.  After  completing  this  process  the  wavefront  is  propagated 
along  the  same  path,  allowing  each  point  to  move  at  a  speed  proportional  to  the  index 
at  the  point  in  the  crystal  in  which  it  finds  itself.  After  a  short  time  the  process  of 
wavefront  movement  is  stopped  and  the  direction  of  each  piece  of  the  wavefront  is 
calculated  using  Huygen's  principle.  The  original  indices  of  the  points  within  any  given 
piece  are  then  propagated  in  that  new  direction,  imparting  those  values  along  the  new 
path  of  the  piece  of  the  wavefront  and  the  process  begins  again.  The  result  of  many 
iterations  of  this  algorithm  is  a  graphic  output  of  the  steady  state  internal  beam 
pattern  due  to  the  photorefractive  effect. 

2.  Limitations 

There  are  several  limitations  to  the  model  that  must  be  discussed  prior  to  the 
introduction  of  results.  The  first,  and  major,  limitation  is  that  of  necessity  each  point 
on  the  wavefront  cannot  be  treated  individually.  The  wavefront  must  be  treated  in 
sections  to  determine  an  appropriate  direction  for  a  section  and  this  makes  the  model 
only  an  approximation  to  the  actual  continuous  case.  The  choice  of  the  number  of 
sections  is  primarily  dependent  on  the  chosen  index  profile.  A  profile  similar  to  that 
shown  in  Figure  3.1  may  be  treated  as  a  single  piece  since  the  direction  of  each 
individual  piece  will  be  the  same  as  the  direction  of  the  entire  wavefront.  By  contrast 
compare  the  profile  of  Figure  3.2.  This  profile  appears  to  lend  itself  to  the  three  pieces 
shovi-n;  however,  in  using  only  three  pieces  the  effect  around  the  extremum  are 
minimized  and  it  may  be  wise  to  consider  the  wave  in  five  pieces. 

The  difference  in  considering  three  as  opposed  to  five  pieces  in  the  wavefront 
may  be  seen  in  Figures  3.3  and  3.4.  In  Figure  3.3  the  index  profile  shown  in  Figure  3.2 
was  used  in  three  pieces.  Note  the  distinctness  of  the  three  discrete  sections.  Figure 
3.4  shows  the  result  of  the  same  profile  divided  into  five  sections.  It  is  still  easy  to 
r.oce  :he  ;Ive  '.ndiviaual  sections  under  consideration,  however  the  output  is  noticeabiy 
different.  Ideally  it  would  be  desirable  to  have  a  very  large  number  of  pieces  used  in 
wavefront  calculations.  However,  the  larger  the  number  of  pieces,  the  larger  the 
number  of  points  in  the  wavefront,  and  consequently  the  larger  the  beam  size.  To  keep 
the  beam/crystal  ratio  in  the  proper  perspective  it  is  necessary  to  increase  the  size  of 
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Figure  3.1  Linear  index  profile  considered  as  one  piece. 

the  crystal  whenever  increasing  the  size  of  the  beam,  and  the  maximum  crystal  size  is 
limited  by  the  available  memoiy'.  The  machine  used  for  this  research  was  limited  to 
two  megabytes  of  random  access  memory,  which  allowed  for  a  maximum  crystal  size  of 
500  units.  To  keep  appropriate  beam/crystal  ratios,  the  beam  size  used  (generally)  was 
between  20  and  40  units  and  the  number  of  pieces  in  the  wavefront  ranged  from  one  to 
eight,  depending  on  the  index  profile. 

A  second  limiting  factor  in  the  use  of  this  simulation  is  seen  in  cases  w’here 
there  is  beam  crossover.  This  happens  when  the  index  profile  has  a  larger  gradient  on 
the  side  away  from  the  c-axis  than  it  does  on  the  c-axis  side.  An  example  of  this  type 
of  profile  can  be  seen  in  Figure  3.3.  In  this  case  the  beam  intensity  profile  is 
-onsiderod  :o  oe  Gaussian  and  the  more  intense  section  of  iie  beam  (i.e.  tiie  portion 
nearest  the  center)  dominates  over  the  weaker  section.  This  ignores  any  consideration 
of  two  beam  coupling  which  w'ould  actually  occur. 

The  final  major  limiting  factor  is  the  inability  of  the  model  to  fully  account  for 
a  piece  of  the  wavefront  which  breaks  away  from  the  main  portion  of  the  incident 
beam.  A  case  w'hcre  this  happens  is  seen  in  Figure  3.3  .  In  this  case  the  index  of 
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Figure  3.2  Nonlinear  profile  with  possible  divisions. 

refraction  gradient  across  the  piece  is  the  same  as  when  it  was  part  of  the  entire  beam. 
This  may  be  a  proper  assumption  if  the  beam  was  continuous  during  fanning;  however, 
as  can  be  seen  in  Figure  2.4  and  many  others  presented  here,  there  are  often  portions 
of  the  beam  which  leave  the  major  portion  of  the  incident  beam  entirely.  Under  the 
condition  of  a  complete  disconnection  it  is  possible  that  the  total  index  gradient  should 
be  based  on  and  not  on  the  index  gradient  as  it  was  calculated  when  'the 

wavefront  piece  was  within  the  incident  beam;  this  was  not  programmed  iiito  the 
model.  There  are  presently  no  data  available  indicating  the  effect  of  a  detached  piece 
of  a  beam  in  contrast  to  the  continuous  case.  This  problem  is  circumvented  when  a 
profile  such  as  that  shown  in  Figure  3.1  is  used  so  that  only  one  wavefront  piece  is 
:ons;dored. 

D.  RESULTS 

As  noted  above,  the  output  of  the  simulation  is  dependent  on  the  index  profile, 
the  beam/crystal  ratio,  the  angle  uf  incidence  and  the  position  on  the  cr>stal  where  the 
incident  beam  enters.  The  criteria  for  success  (of  possible  SPPC)  of  a  given 


Figure  3.3  PSIM  output  for  profile  shown  in  Figure  3.2. 

configuration  was  taken  to  be  the  ability  of  a  portion  of  the  beam  to  reflect  olT  the 
ciystal  corner  in  such  a  way  as  to  rcintersect  itself.  This  would  allow  for  the 
production  of  auxiliary  beams  which  are  a  prerequisite  to  SPPC.  Therefore,  a  result 
such  as  that  seen  in  Figure  3.5  would  be  considered  capable  of  SPPC,  and  one  such  as 
that  seen  in  Figure  3.6  would  not. 

Since  the  purpose  of  this  investigation  is  to  model  beam  fanning,  the  most 
important  parameter  is  the  inde.x  profile.  To  concentrate  on  this  parameter,  a 
beam.'cr\'stai  ratio  and  angle  of  incidence  that  matched  one  of  the  configurations  ^ocn 
in  the  laboratory  were  chosen.  The  configuration  chosen  was  an  often  used 
configuration  with  a  beam/crystal  ratio  of  0.064  and  an  incident  angle  of  40°  .  The 
incident  beam  entered  the  crystal  at  a  point  way  across  the  face  of  the  ciystal.  This 
corresponds  to  a  ciystal  observed  in  the  laboratory  of  5.0  mm  on  a  side  with  a  beam  of 
size  0.32  mm,  incident  3  mm  from  the  front  side  of  the  crystal. 


Figure  3.4  PSIM  output  for  profile  shown  in  Figure  3.2 
but  broken  into  five  separate  pieces. 

Modeling  was  run  with  three  index  profiles.  The  rationale  behind  the  use  of  the 
profiles  will  be  explained  in  Section  E  of  this  chapter.  In  each  case  the  values  of 
^*^max  caused  successful  SPPC  due  to  beam  fanning  are  presented. 

Following  the  theory  of  Feinberg  (Ref.  2],  the  first  index  profile  used  was  that 
shown  in  Figure  2.5.  Successful  SPPC  from  beam  fanning  was  seen  (Figure  3.8), 
however,  note  the  beam  fanning  in  'he  direction  opposite  the  c-axis;  'his  is  never 
observed  in  practice. 

The  second  beam  profile  used  is  shown  in  Figure  3.1.  A  typical  successful  run  is 
shown  in  Figure  3.9.  Although  this  approximates  some  observations  of  beam  fanning, 
note  the  lack  of  a  continued  incident  beam  as  one  often  (but  not  always)  observes. 


figure  3.5  Beam  fanning  resulting  in  successful  SPPC. 

The  final  profile  used  is  showm  in  Figure  3.7.  SPPC  due  to  beam  fanning  was 
exhibited  as  is  seen  in  Figure  3.10.  It  was  this  configuration  which  exhibited  results 
closest  to  those  observed  in  the  laboratory. 

E.  DISCUSSION 

The  first  index  profile  is  taken  directly  from  Feinberg  [Ref.  2].  To  arrive  at  a 
profile  such  as  this  it  is  necessary  to  assume  that  charge  migration  in  the 
photoreiracter  is  symractne.  resulting  .n  a  net  ciiarge  distribution  such  as  siiown  in 
Figure  2.5.  The  beam  pattern  predicted  by  the  simulation  is  shown  in  Figure  3.8  .  The 
major  problem  with  this  profile  is  the  fanning  of  the  beam  in  the  direction  opposite  the 
c-axis  due  to  the  tail  of  the  index  profile.  In  his  original  work  Feinberg  only  considers 
the  central  region  since  it  is  the  most  intense  part  of  the  beam,  ignoring  the  tails  of  the 
profile.  As  can  be  seen  from  the  model,  the  tails  of  the  profile  cannot  be  totally 


Figure  3.6  Beam  fanning  not  resulting  in  successful  SPPC. 


ignored.  To  arrive  at  the  observed  effect  with  the  index  profile  of  Figure  2.5  it  is 
necessaiy  to  assume  that  barium  titanate  exhibits  asymmetric  self-focusing  (i.e.  self- 
focusing  only  in  the  direction  of  the  c-axis).  If  asymmetric  self-focusing  was  present, 
then  the  (less  intense)  light  subject  to  the  profile  of  the  tails  would  be  unable  to  fan 
against  the  c-axis.  To  date  there  has  been  no  reported  evidence  of  self-focusing  of  any 


kind  in  BaTiOji  however,  this  does  not  preclude  the  possibility  since  it  would  only 
affect  a  weak  portion  of  the  beam  and  niav  not  be  obsious  to  b'.e  casual  bi'-'crs er. 
\ote  Wiat  tiiC  I'anning  against  "he  axis  on  :he  .eit-'uand  ude  al'  uie  '^cani  a-.'uiu  'c 
abated  by  two  beam  coupling  and  might  not  be  readily  discerned. 

An  alternative  explanation  is  that  the  crystal  response  to  the  beam  is  not 
characterized  by  the  curves  shown  in  Figure  2.5.  Instead,  all  available  charge  may  he 
totally  excluded  from  all  portions  of  the  beam,  resulting  in  the  profiles  shown  in  Figure 
3.11. 


B«aai  Cantar 

Figure  3.7  Index  profile  for  simulation  of  beam  fanning  in  BaTiO^. 

In  an  effort  to  approximate  an  index  profile  which  would  exhibit  the  necqssarx' 
behavior,  the  profile  shown  in  Figure  3.1  was  used.  This  beam  profile  was  arrived  at 
by  making  the  important  assumption  that  there  is  a  total  exclusion  of  charge  carriers 
within  the  illuminated  section  of  the  cr>'stal.  This  is  an  easy  profile  to  work  with  sirree 
the  outcome  is  independent  of  the  number  of  pieces  considered  in  the  wavefront,  as 
previously  mentioned.  The  result  of  the  simulation  with  this  profile  is  shown  in  Figure 
3.9.  However,  there  are  two  major  problems  with  this  profile.  The  first  one  is  the  lack 
of  extension  of  the  incident  beam  as  mentioned  above.  Note  that  in  most  of  the 
photographs  presented  in  Appendix  A  there  is  an  extension  of  the  incident  beam.  This 
is  almost  universally,  though  not  exclusively,  true  (see  Figure  .A. 5  for  c\a tuple f 
•Second,  it  is  unrealistic  ;n  its  siiarp  -ransitions.  ihu;;  )i'-;icse  prci'icn'.s  .ro  •  ii'.cu 
the  introduction  of  the  third  inucx  proiiie. 

The  final  index  profile  (Figure  3.7)  resulted  in  a  beam  geometry  which  exhibited 
all  of  the  requisite  properties.  This  profile  was  arrived  at  by  the  rounding  of  the 
discontinuities  of  the  prior  profile;  specifically,  a  sine  function  was  used  to  arrive  at  the 
configuration  shown  in  Figure  3.7.  The  results  of  the  siinulation  are  presented  in 


Figure  3  i  PSIM  output  using  the  index  profile  of  Figure  2.5 

Figure  3  10  To  arrr.e  at  this  profile  it  '.vouid  be  necessary  for  difiusion  of  the  charge 
carriers  almost  completely  out  oi  the  incident  beam  path.  The  eventual  dilfusion  of 
charge  carriers  completely  out  of  the  path  may  actually  be  expected  since  anv  light  of 
intensity  above  ambient  'mH  cause  any  charge  curriers  to  eventually  be  totally 
liberated'*  Die  graphs  of  Figure  5  11  show  the  ncccs'-ary  charge,  electric  I'leld  and 
intensity  profile  ot  the  incident  beam  necessary  to  arrive  at  an  index  profile  similar  to 


appearance  of  tl.e  c.mitinucJ  beam. 


■*1  einbcrc  stiows  hov.evcr.  tliat  the  electric  field  may  be  proportional  to  tl.e 
gradient  cl  the  in..Kle;it  bc.im  mien',;".  'Ref  2'! 
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To  arrive  at  the  configuration  of  Figure  3.9  a  value  of  An^^ay  *  0.03  was  used. 
In  Figure  3.10  -  0.02.  Regardless  of  the  profile  used,  these  values  of  An^j^^^^ 

are  characteristic  of  the  value  necessary  to  account  for  the  observed  beam  fanning  in 
BaTiOj.  For  successful  SPPC  with  the  configuration  chosen,  An^j^j^^  must  be  greater 
than  10'^.  Reducing  the  beam  size  will  obviously  reduce  the  necessary  value  of 
Aumax-  values  bn  the  order  of  10*^  to  10*^  are  typical. 

Thus,  the  simulation  has  predicted  the  necessary  change  in  index  of  refraction 
across  the  beam  path,  and  the  necessary  profile  of  that  change.  In  Chaper  V  the 
results  of  the  necessary  value  for  will  be  compared  to  those  possible  for  the 

cry'stal.  Also,  the  physical  implications  of  the  necessary  profile  of  the  index  of 
refraction  within  the  beam  path,  as  determined  by  the  simulation,  will  be  discussed. 
The  comparison  of  the  information  derived  from  the  simulation  to  that  experimentally 
denved  from  the  crystal,  will  show  that  the  necessary  conditions  cannot  be  met  in 
BaTiOj. 


IV.  EXPERIMENTS 


A.  INTRODUCTION 

In  order  to  investigate  beam  fanning  in  BaTiO^,  an  experiment  was  designed  to 
measure  the  influence  of  a  change  in  the  incident  beam  diameter  on  the  fanning 
process.  If  fanning  increases  with  a  decrease  in  beam  diameter,  it  could  be  concluded 
that  beam  farming  is  a  direct  result  of  the  change  in  index  of  refraction  within  the 
beam  path,  as  predicted  in  Equation  2.25.  If,  however,  the  fanning  decreases  with  a 
decrease  in  incident  beam  diameter,  it  could  be  concluded  that  beam  fanning  is  actually 
forward- scattering  due  (probably)  to  SPS,  since,  as  shown  in  Equation  2.28,  the  gain 
per  unit  length  will  decrease  with  a  decrease  in  spot  size. 

B.  EXPERIMENTAL  APPARATUS 

The  apparatus  used  for  this  work  is  shown  in  Figure  4.1.  The  laser  was  a  Spectra 
Physics  Ar"*"  variable  power  (20  mW  maximum),  multimode  laser,  model  #162A-07. 
During  the  experimentation  the  laser  was  operated  at  all  of  the  visible  lines  and  was 
not  optically  isolated  from  the  experiment.  Upon  leaving  the  laser,  the  beam  was 
directed  into  a  polarization  rotator  in  order  to  rotate  the  beam  from  the  ordinary  to 
the  extraordinary  polarization  in  the  crystal.  To  ensure  that  only  extraordinarily 
polarized  light  was  used,  the  beam  was  then  directed  through  a  Glan-Thompson 
polarizer.  Following  the  polarizer  was  a  beam  splitter,  to  enable  the  observation  of 
any  phase-conjugate  reflection,  and  a  variable  aperture  of  the  type  common  in 
photographic  lenses.  This  aperture  acted  as  a  spatial  filter,  reducing  the  amount  of 
scattered  light  reaching  the  cr>stal.  The  beam  was  then  focused  onto  the  crystal  by 
uncoaied  lenses  of  either  34  cm  or  16  cm  focal  length. 

The  crystal  was  one  of  four  single  crystals  of  BaTiOj  manufactured  by  Sanders 
•Associates  and  poled  into  a  single  domain.  Two  of  the  crystals  ion  loan  from  the  Los 
.\iamos  National  Laboratories)  measured  approximately  5  x  .5  x  5  mm  and  were 
transparent  with  a  slight  yellow  tint.  Two  of  the  crystals  (belonging  to  the  Naval 
Postgraduate  School)  measured  approximately  5  x  5  x  2.5  mm  and  had  no  noticable 
tint  to  them.  The  crystal  was  mounted  on  a  rotatable  platform  capable  of  being 
translated  12  cm  along  the  axis  of  the  optical  bench  on  which  all  of  the  components 
were  mounted.  Following  the  crystal,  and  connected  to  the  mounting  table,  was  a 


knife  edge.  The  knife  edge  was  coated  with  black  tape  to  ensure  that  a  resonating 
cavity  was  not  formed  with  the  crystal.  Behind  and  to  the  side  of  the  knife  edge,  also 
connected  to  the  translation  table,  was  a  silicon  photodiode  connected  to  an 
electrometer  measuring  the  short  circuit  current. 

C.  EXPERIMENTAL  PROCEDURE 

To  measure  the  amount  of  beam  farming  quantitatively,  the  knife  edge  was 
aligned  in  such  a  way  that  the  beam  transmitted  through  the  crystal  was  blocked  from 
the  detector.  To  ensure  that  the  blockage  of  the  transmitted  beam  was  sufficient,  the 
table  w'as  translated  rapidly  along  its  entire  range  of  motion  and  the  detector  output 
was  monitored  to  ensure  that  there  was  little  or  no  change  in  the  recorded  intensity. 
The  vibration  inherent  in  the  movement  and  the  changing  beam  diameter  of  the 
incident  beam  ensured  the  lack  of  significant  beam  fanning,  so  that  a  lack  of  change  in 
the  recorded  intensity  indicated  that  the  transmitted  beam  was  completely  blocked 
from  the  detector  for  all  observed  beam  diameters. 

Once  the  ciy'stal,  knife  edge  and  detector  were  in  position,  the  beam  was  blocked 
and  the  crystal  was  flooded  with  light  from  a  75  watt  incandescent  lamp  located 
approximately  12  cm  away.  After  approximately  one  minute  the  light  was  turned  off 
and  the  laser  beam  was  unblocked.  The  intensity  of  the  light  reaching  the  detector  was 
then  monitored  and  the  peak  and  steady-state  intensities  were  recorded.  Once  the 
intensity  had  reached  a  steady-state  value  the  beam  was  again  blocked,  the  crystal  was 
translated  and  flooded  with  white  light,  and  the  procedure  was  repeated. 

The  flooding  of  the  crystal  with  white  light  effectively  homogenized  the  charge 
distribution  within  the  crystal  and  ensured  that  one  reading  was  not  affected  by  the 
previous  one.  Special  care  w’as  taken  to  ensure  that  the  beam  was  parallel  to  the 
crystal  table  and  that  all  components  were  mounted  on  the  optical  axis  of  the  optical 
bench.  The  precision  translation  mechanism  on  which  the  crystal  table  was  mounted 
ensured  that  the  beam  was  incident  on  the  crystal  with  the  same  center  for  all 
measurements.  These  precautions  effectively  compensated  for  any  surface 
-.nconsistcnc-.es  .on  '.he  cr.'stat  face  except  at  ver/  .small  laser  beam  diameters. 

Both  peak  and  steady-state  values  of  the  fanned  beam  intensity  were  recorded; 
however,  only  the  peak  values  were  used  for  final  analysis.  This  decision  was  made 
because,  although  the  steady  state  values  generally  correlated  with  the  peak  values,  the 
steady  state  values  could  be  affected  by  numerous  variables.  One  obvious  problem 
with  steady-state  readings  is  the  possibility  of  phase  conjugation,  thus  reducing  the 


amount  of  fanning  recorded  in  the  steady  state.  Also,  due  to  the  orientation  of  the 
detector,  it  was  possible  (even  likely)  that  the  beam  would  fan  across  the  detector, 
eventually  coming  to  rest  past  the  area  it  can  record;  sometimes  the  fanning  was 
observed  to  exit  the  face  perpendicular  to  the  face  the  detector  was  monitoring.  To 
reduce  all  of  these  effects,  the  incident  beam  was  originally  allowed  to  enter  at  only 
small  angles  to  the  c-axis  (*^2®  -  5®)  and,  although  alignment  was  tedious,  consistent 
results  were  obtained  for  various  crystals,  lens  focal  lengths,  laser  lines  and  incident 
angles.  Eventually  larger  incident  angles  were  used  (-^20®  -  40®)  with  similarly 
consistent  results. 

To  ensure  that  the  observed  effect  is  not  influenced  by  scattering  initiated  by 
small  (e.g.  molecular  size)  crystal  inhomogeneities,  one  configuration  was  chosen  in 
which  there  was  moderate  fanning  and  the  fanned  intensity  was  recorded  for  each  of 
the  available  visible  lines  of  the  laser.  The  result  was  a  mean  ratio  of  incident  intensity 
to  fanned  intensity  that  was  constant  across  all  available  wavelengths,  thus  indicating 
no  resonances  and  no  dependence  of  scattering  on  wavelength  at  the  available  laser 
lines. 

After  the  data  were  collected,  the  beam  diameter  at  each  point  was  measured 
with  an  Aeronca  Electronics  Laser  Blade  and  the  direct  laser  intensity  was  recorded. 
Here,  the  beam  diameter  is  defined  as  twice  the  distance  from  the  center  of  the  beam  to 
the  point  of  10%  intensity.  For  analysis  purposes  the  beam  diameter  could  be  divided 
by  the  cosine  of  the  angle  of  incidence  so  that  the  beam  diameter  would  be  that 
actually  along  the  y-axis  of  the  crystal. 

Measurements  were  also  made  by  placing  the  detector  above  the  crystal  as  the 
incident  beam  diameter  was  changed.  This  gave  an  indication  of  the  total  amount  of 
scattering,  not  just  asymmetric  scattering.  In  this  case,  the  front  surface  of  the  crystal 
was  shielded  from  the  detector  with  black  tape  to  keep  specular  reflection  off  of  surface 
irregularities  from  influencing  the  measurements.  Transmission  measurements  were 
also  made  by  replacing  the  knife  edge  with  an  aperture,  thus  allowing  only  the 
transmitted  portion  of  the  incident  beam  to  reach  the  detector. 

D.  RESULTS 

Typical  results  of  the  fanning  measurements  are  shown  in  Figure  4.2.  In  this 
case  the  crystal  was  at  an  angle  of  <  5®  so  the  beam  diameter  need  not  be  adjusted  for 
COS0.  All  intensity  recordings  have  been  normalized  for  comparison  purposes.  The 
shape  of  the  curves  is  insensitive  to  the  laser  line  used,  beam  angle,  choice  of  crystal, 
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and  the  divergence  or  convergence  of  the  incident  beam  (i.e.  whether  the  crystal  is 
placed  in  front  of  or  behind  the  focal  point  of  the  focusing  lens). 

Figure  4.3  is  a  plot  of  the  intensity  seen  from  above  the  crystal  for  various  beam 
diameters  for  converging  light.  There  are  two  important  things  to  note  here.  First, 
there  is  less  scattering  for  smaller  beam  diameters  and  second,  the  sharp  discontinuity 
at  about  three  micron  beam  diameter.  This  discontinuity  is  noticable  as  the  crystal 
goes  through  the  focal  point  of  the  focusing  lens.  Figure  4.4  is  a  plot  of  the  scattering 
as  viewed  from  above  the  crystal  for  diverging  light.  Note  the  initial  decrease  in  the 
scattering  within  the  first  four  microns.  This  is  only  apparent  with  diverging  light  and 
is  seen  in  all  crystals  I  used.  Note  also  that  there  is  a  much  broader  dynamic  range  in 
the  intensity  for  diverging  light  than  for  converging  light. 

Although  all  of  the  scattering  data  are  generally  insensitive  to  which  crystal  is 
used  (i.e.  the  general  shape  of  all  the  curves  for  a  given  situation  is  the  same  for  all 
crystals),  the  same  cannot  be  said  for  the  transmitted  beam.  Figures  4.5  and  4.6  show 
the  transmitted  intensity  for  all  four  crystals  plotted  together.  Crystal  #2  is  the  only 
one  which  does  not  exhibit  SPPC  and  is  the  only  crystal  in  which  the  transmitted  beam 
intensity  increases  with  beam  diameter  as  a  general  trend.  However,  note  that  all  of 
the  curves  exhibit  some  inflection  point  around  six  micron  diameter  for  diverging  light. 
This  inflection  is  not  a  consistent  maximum  or  minimum;  it  depends  on  the  crystal. 
Presently  I  cannot  explain  the  cause  of  these  dissimilar  curves.  Generally  it  may  be 
noted  that  the  transmitted  intensities  map  the  same  characteristics  as  the  scattered 
intensity,  but  inverted,  as  would  be  expected.  All  of  these  various  features  will  be 
addressed  later,  but  the  important  aspect  of  all  of  the  data  is  obvious;  increased  beam 
diameter  results  in  increased  beam  fanning. 

Although  these  data  alone  are  a  strong  indication  of  the  lack  of  index  change 
w'ithin  the  beam  path  due  to  the  photorefractive  effect,  there  exists  even  stronger 
evidence.  By  removing  the  knife  edge  and  detector  from  behind  the  crystal,  the  image 
of  the  transmitted  spot  could  be  observed  on  a  screen  approximately  three  meters 
lu-a-.  <ho’.vs  -he  image  of  the  spot,  which  was  projecied  directly  onto 

pnctograpmc  .:im  a'  aiir'ercnt  times  after  the  beam  was  allowed  to  enter  the  cr^/stai, 
and  prior  to  the  onset  of  phase  conjugation.  Figure  4.7a  is  an  image  of  the  beam 
immediately  after  turn  on;  notice  the  spot  is  generally  undistorted.  Figure  4.7b  is  a 
photograph  between  turn  on  and  steady-state,  and  Figure  4.7c  is  a  photograph  at 
steady-state.  By  viewing  these  images  it  is  clear  that  the  energy  of  the  beam  is  first 
extracted  from  the  outside  of  the  beam  and  not  from  the  inside. 
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In  addition  to  the  above  measurements,  the  index  of  refraction  of  the  crystal  was 
estimated  by  measuring  Brewester's  angle.  At  X  *  514  nm  the  published  value  for  the 
extraordinary  index  of  refraction  is  n^  =  2.424  [Ref  22],  however,  my  measurements 
consistently  resulted  in  n^  =  2.26  ±0.02.  This  large  discrepancy  is  still  unexplained.  It 
is  also  interesting  to  note  that,  although  there  was  a  distinct  minimum  at  Brewester's 
angle,  the  reflection  did  not  disappear  entirely.  The  minimum  was  so  slight  that  it  was 
quite  difficult  to  find,  and  eventually  a  photodetector  had  to  be  used  to  accurately 
locate  it.  It  is  possible  that  both  of  these  observations  may  be  explained  by  the 
presence  of  impurities  on  the  surface  of  the  crystal.  Specifically,  an  oil  film  from  the 
hands  of  the  experimenter  may  be  at  fault.  Additionally,  there  is  the  possibility  that 
heating  at  the  surface  of  the  crystal  due  to  the  high  energy  density  of  the  incident  beam 
may  be  responsible.  Since  measurements  by  others  [Ref  30]  of  the  index  of  refraction 
of  BaTiOj  using  Brewester's  angle  have  resulted  in  indices  consistent  with  the  reported 
values,  either  one  or  both  of  these  explanations  may  account  for  the  observed 
anomalies. 
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Figure  4.7a  Photograph  of  projected  spot  immediately  upon  turning  on  the  incident  beam. 


V.  ANALYSIS 


A.  DISCUSSION  OF  RESULTS 

The  simulation  in  Chapter  III  and  the  experiments  described  in  Chapter  IV 
clearly  show  that  asymmetric  self-defocusing  cannot  be  due  solely  to  a  change  in  index 
of  refraction  along  the  beam  path  within  the  crystal.  Even  with  the  changes  of  the 
profile  for  the  change  in  index  of  refraction  proposed  in  Chapter  III,  the  necessary 
value  of  ^^8®  according  to  the  published  maximum  values  for  BaTiOj. 

Inspection  of  the  necessary  values  of  An^^^  from  the  simulation  shows  that  values  on 
the  order  of  10*^  are  necessary  to  account  for  the  observed  fanning  within  the  crystal; 
however,  previously  calculated  values  place  an  upper  bound  of  ^  10*^,  and  a  value  of 
~  10"^  is  probably  more  reasonable^  (Ref.  11]. 

This  discrepancy  may  be  accounted  for  only  by  assuming  that  the  index  gradient 
is  limited  to  a  very  small  portion  (in  some  cases  as  small  as  '/jq  the  size  of  the  input 
beam)  of  the  center  of  the  beam,  or  by  assuming  that  the  observed  effea  is  a  result  of 
another  process.  It  is  impossible  to  justify  the  first  assumption  and  still  achieve  the 
necessary  index  profile,  since  charge  carriers  must  be  liberated  from  the  entire  region 
where  the  brightness  of  the  incident  beam  exceeds  ambient.  Therefore,  by  Laplace's 
equation,  the  gradient  of  the  electric  field  within  the  beam  path  must  be  uniform.  This 
leaves  only  the  possibility  that  another  process  is  responsible  for  beam  fanning. 

The  results  of  the  experiments  of  Chapter  IV  reinforce  the  conclusions  made 
from  the  simulation  portion.  Figure  4.2  shows  that  beam  fanning  decreases  with  a 
decrease  in  beam  diameter.  From  Equation  2.25  it  is  clear  that  this  behavior  cannot  be 
explamed  by  a  change  in  the  index  of  refraction  within  the  beam  path  due  to  the 
photorefractive  effect.  If  this  type  of  change  were  responsible  for  beam  fanning,  the 
the  amount  of  fanning  would  have  to  increase  with  a  decrease  in  beam  diameter. 
X.ivi;:; :r.ail*. .  howr.  n  7 '’urc  'eram  v.nninc  •jeems  '^n  ‘he  "'or.rjher.’  lad  'he 
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the  index  profile  as  shown  in  Figure  2.2,  this  is  impossible  if  the  fanning  is  due  to  an 
asymmetric  change  in  index  of  refraction  across  the  beam  path,  as  was  originally 
predicted. 

^Feinberg  does,  however,  calculate  a  possible  maximum  value  on  the  order  of 
10*^  (Ref  2^]. 


One  may  argue  that  all  of  the  observed  qualities  of  beam  fanning  may  be 
explained  by  assuming  "small-scale"  self-focusing  of  the  incident  beam  within  the 
crystal.  That  is,  the  incident  beam  may  experience  a  breakup  into  spots  which  are  of 
diameters  small  enough  to  create  index  gradients  of  the  magnitude  shown  necessary  by 
the  simulation  in  Chapter  III.  This  type  of  self-focusing  has  been  known  for  over  two 
decades,  and  is  common  in  materials  whose  dielectric  constant  increases  with  field 
intensity  (Ref  31].  If  small-scale  self-focusing  occurred  in  BaTiOj,  it  could  explain 
many  of  the  effects  that  are  observed  during  SPPC  (e.g.  the  presence  of  small  filaments, 
beam  fanning,  etc.).  However,  there  are  serious  problems  with  a  small-scale  self- 
focusing  theory.  First,  the  breakup  of  optical  beams  due  to  self-focusing  has  been 
shown  to  be  spatially  periodic,  with  the  periodicity  being  dependent  upon  intensity 
[Ref  32).  As  is  demonstrated  in  Figure  4.7,  there  is  no  indication  of  spatial  periodicity 
in  the  beam  as  it  is  transmitted  through  a  BaTiO^  crystal.  Also,  even  in  good  self- 
focusing  mediums  at  very  high  beam  intensities  ('^10^  W/cm^),  it  is  necessary  for  the 
beam  to  have  a  long  path  length  (on  the  order  of  tens  of  centimeters)  within  the 
medium  before  self-focusing  is  observed  [Ref  33).  Even  with  a  very  large  change  in  the 
dielectric  constant,  it  is  unlikely  that  self-focusing  could  occur  within  a  medium  of  only 
5  mm  in  length.  Finally,  there  is  no  evidence  in  the  literature  that  the  dielectric 
constant  of  BaTiOj  changes  significantly  with  the  intensity  of  the  incident  light 
(although  it  does  change  with  an  intensity  gradient,  as  discussed  in  Chapter  II).  It  is 
interesting  to  note  however,  that  the  dielectric  constant  of  BaTiO^  does  change  with 
temperature  [Ref  25],  and  this  may  cause  self-focusing  at  very  high  powers  and  very 
small  beam  diameters,  given  a  long  enough  path  within  the  medium. 

From  the  evidence  presented  here  it  follows  that  asymmetric  self-defocusing  in 
BaTiOj  cannot  be  due  to  an  asymmetric  change  in  index  of  refraction  along  the  path 
of  the  incident  beam  due  to  the  photorefractive  effect.  Attributing  this  phenomenon  to 
some  other  process  however,  is  not  so  easy. 


B.  SC.kTTERING 
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explain  asymmetric  self-defocusing  in  BaTiOj,  besides  the  theory  just  discussed,  is  the 
theory  involving  scattering  (SPS)  [Ref  24].  Although  I  have  not  conclusively 
demonstrated  that  beam  fanning  in  BaTi03  is  due  to  SPS,  there  are  a  number  of 
observations  which  indicate  that  this  type  of  scattering  maybe  the  origin  of  the  effect. 
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As  may  be  seen  in  Equation  2.28,  if  near-forward  scattering  is  the  origin  of  beam 
fanning,  then  there  should  be  a  decrease  in  fanning  with  a  decrease  in  incident  beam 
diameter,  an  observation  borne  out  in  our  experiments.  Also,  the  diversion  of  energy 
from  the  incident  beam  could  progress  from  the  periphery  toward  the  center,  since 
noise  photons  originating  from  outside  the  beam  must  interact  with  the  first  portion  of 
the  beam  they  encounter  (Figure  4.7). 

Another  interesting  and  important  discovery  is  that  although  the  shape  of  the 
fanned  intensity  curves  is  similar  regardless  of  convergence  or  divergence  of  the 
incident  beam  (in  front  of  or  behind  the  focal  point  of  the  focusing  lens),  the  absolute 
intensity  of  the  scattered  light  is  always  observed  to  be  greater  for  the  diverging  case. 
This  is  not  perplexing  when  considering  near-forward  scattering.  One  would  expect 
more  scattering  from  diverging  light  than  converging  because  the  size  of  the  beam 
increases  as  it  progresses  through  the  crystal  for  diverging  light,  thus  increasing  the 
amount  of  scattering.  Also  it  is  noted  that  individual  intensity  readings  are  not 
consistent,  even  if  there  is  no  no  movement  of  the  apparatus  between  readings.  This  is 
expected  since  the  scattering  process  begins  with  noise  photons,  and  no  two  situations 
can  be  exactly  alike,  even  for  similar  configurations.  This  inability  to  exactly  reproduce 
any  given  measurement  suggests  that  the  scattering  may  begin  with  random  index 
changes  within  the  crystal,  caused  by  the  stochastic  distribution  of  charge  carriers. 

Finally,  if  stimulated  scattering  is  responsible  for  beam  fanning,  then  one  would 
not  expect  continuous  fanning,  but  fanning  only  along  paths  from  which  an  incident 
photon  interacted  with  the  incident  beam,  thus  creating  a  set  of  filaments  emanating 
from  the  incident  beam.  A  close  examination  of  any  of  the  photographs  of  Appendix 
A  show’s  that  this  is  indeed  the  case.  Fanning  is  a  set  of  discrete  small  beams  all 
emanating  from  the  incident  beam.  We  must  now  discuss  curved  appearance  of  the 
dominant-scattered  beam  during  self-pumped  phase  conjugation,  a  phenomenon  not 
explained  by  SPS. 

A  very  close  examination  of  many  photographs  of  SPPC  (e.g.  Figure  A. 5)  reveals 
'hat  beams  ’vhich  appear  to  cur.’e  often  have  a  series  of  discontinuities  within  '.hem. 
Generaily,  a  rilament.  or  set  of  filaments,  leaves  '.he  incident  beam  along  'he  direction 
of  the  c-axis.  After  a  certain  distance  (usually  about  halfway  to  the  crystal  edge)  these 
filaments  enter  a  region  where  they  become  fuzzy,  and  often  two  become 
indistinguishable.  The  filaments  exit  this  fuzzy  region  at  an  angle  different  than  their 
entrance  angle.  This  type  of  behavior  is  explainable  by  assuming  that  within  this  fuzzy 


region  there  is  stimulated  scattering  due  to  SPS.  Forward  scattered  beams  enter  into  a 
region  where  scattering  occurs  due  to  photons  arriving  from  a  different  angle.  The 
beams  which  continue  on  at  a  new  angle  are  the  forward  scattered  beams  from  this 
second  scattering.  It  is  therefore  a  series  of  this  type  of  scattering  discontinuities  which 
gives  the  characteristic  curt'ed  appearance  of  a  filament  during  SPPC.  However,  this  is 
not  always  the  case. 

There  exist  configurations  where  this  series  of  discontinuities  does  not  appear  to 
account  for  the  curved  filaments  (i.e.  Figure  A.6).  In  a  majority  of  these  cases  a  close 
examination  reveals  that  a  large  filament  originates  at  the  incident  beam  or  a  series  of 
filaments  leave  the  incident  beam  very  close  together.  The  curved  appearance  comes 
from  the  coupling  of  energy  to  the  filament,  or  portion  of  the  filament,  in  the  direction 
of  the  c-axis  as  the  filaments  progress  through  the  crystal.  Therefore  the  beam  doesn't 
actually  curve;  it  simply  appears  to  do  so  because  it  is  losing  energy  from  one  side  and 
transferring  it  to  the  other.  The  process  is  actually  a  form  of  self-focusing  due  to  two 
beam  coupling. 

I  have  observed  a  veiy’  few  instances  where  it  appears  that  neither  of  these 
processes  can  fully  account  for  the  curved  appearance  of  the  fanned  filaments.  In 
these  cases  one  would  expect  the  effect  to  be  due  to  true  beam  fanning  from  the 
change  in  index  of  refraction  within  the  filament  path.  Although  the  simulation  PSIM 
shows  that  it  is  impossible  for  this  phenomenon  to  fully  account  for  the  observed 
fanning  effect,  it  would  be  surprising  not  to  find  some  evidence  of  a  change  in  index  of 
refraction  on  the  order  of  10*^  across  a  beam  of  diameter  less  than  a  few  microns. 

C.  SELF  PUMPED  PHASE  CONJUGATION 

From  the  discussion  of  the  previous  section  we  may  conclude  that  a  change  in 
index  of  refraction  within  the  incident  beam  path  is  not  responsible  for  asymmetric  self 
defocusing  in  BaTiOj.  It  is  probable  that  SPS  is  responsible  for  the  effect.  Therefore 
one  may  say  that  the  true  origin  of  SPPC  in  BaTiOj  is  scattering,  probably  SPS.  The 
cause  of  the  final  phase  conjugate  reflection  is  not  as  clear. 

'’base  coniugation  by  stimulated  scattering  in  BafiO,  iias  oeen  demonstrated 
[Ref  23],  and  recently  it  has  been  shown  how  backscattering  may  dominate  forward 
scattering  in  SPS  [Ref  24].  Additionally,  the  investigations  reported  here  indicate  that 
scattering  is  responsible  for  the  phenomenon  of  beam  fanning.  Therefore  one  might 
incorrectly  conclude  that  phase  conjugation  in  BaTiOj  is  due  to  stimulated  scattering. 
With  the  onset  of  backward  scattering  in  the  crystal  there  begins  a  competition  of  the 


scattered  waves.  As  the  modes  begin  to  build,  a  mode  showing  high  gain  can  usurp 
power  from  one  showing  less  gain.  In  this  way  only  the  modes  showing  the  highest 
gain  eventually  survive.  It  is  this  competition  between  configurations  having  high  gain 
which  is  responsible  for  the  amplitude  oscillations  often  seen  in  a  phase  conjugate 
reflections  (see  [Ref.  34,35]  or  Appendix  A). 

If  a  filament,  arising  from  near  forward  scattering,  has  a  configuration  such  that 
it  can  reflect  off  the  comer  of  the  crystal  (i.e.  a  CAT  comer)  to  create  a  situation 
w’here  degenerate  four  wave  mixing  (DFWM)  is  possible  in  addition  to  SPS,  then  one 
would  expect  that  filament  to  have  a  gain  advantage.  This  then  would  be  the  filament 
to  eventually  win  the  competition.  So  it  actually  would  be  expected  that  both  DFWM 
and  SPS  should  contribute  to  SPPC. 

A  review  of  the  conditions  of  the  only  published  report  of  phase  conjugation  due 
solely  to  SPS  [Ref  23]  shows  that  special  precautions  were  taken  to  ensure  maximum 
coupling.  Specifically,  an  index  matching  liquid  was  used  which  preserved  as  much  of 
the  phase  conjugate  return  as  possible.  Additionally  the  crystal  was  surrounded  by  a 
glass  cuvette  (containing  the  index  matching  liquid)  which  increased  the  amount  of 
scattered  light  available  for  SPS  to  begin.  Even  with  these  advantages,  a  phase 
conjugate  reflection  of  only  ~  10%  was  obtained.  This  is  quite  small  compared  to  the 
30%  to  50%  usually  seen  from  SPPC  in  BaTiOj  in  air.  This  indicates  that  DFWM 
from  the  internal  comer  reflection  is  the  most  significant  contributor  to  the  phase 
conjugate  return  observed  in  SPPC,  as  was  originally  proposed  [Ref  2]. 


VI.  CONCLUSIONS 


Experimental  evidence  does  not  support  the  theory  that  asymmetric  self 
defocusing  in  BaTiOj  is  a  result  of  an  asymmetric  change  in  the  index  of  refraction 
within  the  path  of  the  incident  beam.  Instead,  forward  stimulated  scattering,  probably 
Stimulated  Photorefractive  Scattering  (SPS),  is  found  to  be  a  more  likely  explanation. 

Simulation  has  shown  that  for  an  index  change  to  be  responsible  for  the 
phenomenon  of  beam  fanning,  the  necessary  change  in  the  index  of  refraction  across 
the  beam  must  exceed  the  theoretical  maximum,  based  on  the  maximum  possible  space 
charge  field.  A  change  in  the  index  of  refraction  on  the  order  of  10"^  is  necessary  to 
account  for  the  observed  effects,  whereas  the  accepted  maximum  is  on  the  order  of  10“^ 
to  10"^. 

Experimental  results  do  not  support  a  theory  which  requires  a  change  in  index  of 
refraction  across  the  incident  beam  as  the  cause  of  beam  fanning.  The  magnitude  of 
fanning  is  observed  to  increase  with  an  increase  in  beam  diameter.  Also,  observations 
of  the  maimer  in  which  the  beam  fans  show  that  energy  is  scattered  out  of  the  beam 
path  beginning  with  the  periphery  and  not  beginning  in  the  center.  Neither  of  these 
observations  is  consistent  with  the  results  expected  if  a  purely  photorefractive  effect 
were  responsible  for  asymmetric  self  defocusing. 

Close  observations  of  the  process  of  beam  fanning  in  BaTiOj  indicate  that  beam 
fanning  can  be  explained  by  a  stimulated  scattering  theory.  The  presence  of  filaments 
rather  than  continuous  fanning,  scattering  beginning  with  the  periphery  of  the  incident 
beam,  and  discontinuous  filaments  during  phase  conjugation  are  all  explicable  when 
beam  fanning  is  attributed  to  stimulated  scattering. 

Future  experiments  in  this  area  should  include  probing  the  "fuzzy"  regions  of  the 
observed  filaments  from  above  with  a  low  power  laser  to  search  for  the  presence  of 
additional  Inde.'c  gratings.  Also,  the  actual  mange  in  the  index  of  refraction  across  tne 
incident  oeam  path  witiim  the  crv'stai  remains  to  be  determined  experimentally. 
Additionally,  further  research  into  the  size,  type  and  spacing  of  scattering  centers 
within  the  crystal  is  necessary  to  fully  explain  the  observed  difference  in  scattering  of 
diverging  and  converging  light. 
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APPENDIX  A 

INTERNAL  BEAM  PRODUCTION 
1.  INTRODUCTION 

The  production  of  internal  beams  in  BaTiOj  has,  to  my  knowledge,  never  been 
systematically  investigated.  A  search  of  the  available  literature  indicates  that  the  cause 
for  this  oversight  may  be  that,  in  general,  millimeter  sized,  collimated  beams  have  been 
used.  In  all  available  photographs  of  BaTiO^  used  as  a  SPPCM,  the  beam  has  not 
been  focused  onto  the  crystal.  If  the  beam  is  not  focused,  then  the  diameter  of  the 
beam  upon  entering  the  crystal  is  at  best  the  same  as  when  it  leaves  the  laser.  It 
appears  that  when  a  beam  is  not  focused  onto  the  crystal,  auxihary  beams  are 
eliminated  or  obscured  by  the  strong  input  beam  since  its  diameter  is  large  compared 
to  the  dimensions  of  the  crystal  (sometimes  being  more  than  half  the  size  of  the 
crystal).  However,  if  the  input  beam  is  focused,  the  beam  diameter  is  small  compared 
to  the  crystal  and  the  immediate  effects  of  the  input  beam  stay  localized  within  the 
crystal,  thus  allowing  the  formation  and  observation  of  auxiliary  beams.  There  is  also 
an  immense  increase  in  the  energy  density  of  the  beam  within  the  crystal  upon 
focusing. 

It  was  using  a  focused  beam  that  we  first  discovered  the  ability  of  barium 
titanate  to  support  internal  auxiliary  beam  patterns.  Extensive  investigation  by  A.  V. 
Nowak  and  myself  at  Los  Alamos  National  Laboratories  indicated  that  there  is  a 
plethora  of  internal  beams  which  may  be  generated  by  the  nonlinear  effects  of  barium 
titanate  during  phase  conjugation.  It  appears  that  phase  conjugation  is  necessarv’  to 
sustain  these  beams  since  no  au.xiliary  beams  were  observed  in  the  absence  of  phase 
conjugation.  In  many  cases  it  seems  that  the  beam  is  initiated  by  the  reflection  of  the 
phase-conjugate  beam  off  of  the  face  of  the  crystal  as  it  exits,  but  this  is  not  always 
true. 

'•Vhat  follows  :a  ‘.his  chapter  ’.s  an  attempt  to  :-eporL.  .'ategori.’o.  ar.J  :\o!ain 
many  or't.ne  ooser'-'ed  ;nternai  beam  patterns  ana  anomalous  erfects.  It  :s  ;;r.por*..r'.:  'o 

note  that  often  it  is  difficult  or  impossible  to  recreate  some  of  the  observed  effects.  It 
appears  that  not  only  position  and  beam  profile  are  important,  but  also  the  short  term 
history  of  the  crystal.  This  means  that  possibly  one  effect  may  not  be  observable 
unless  another  has  (or  has  not)  immediately  preceded  it. 


2.  EXPERIMENTAL  APPARATUS 


All  investigations  presented  in  this  chapter  were  performed  at  Los  Alamos 
National  Laboratorv’  on  an  air  floated  optical  bench  in  free  air  that  was  maintained  at 
a  temperature  of  22±  1°C.  The  apparatus  used  is  shown  schematically  in  Figure  A.l. 

During  all  experimentation  the  light  source  was  a  Spectra  Physics  2020-03  argon 
ion  laser,  prism-tuned  to  514.5  nm  and  operated  in  a  single  longitudinal  mode  by 
means  of  an  intercavity  eialon.  The  beam  at  the  laser  had  a  diameter  of  0.88  mm  at 
the  of  peak  intensity. 

After  leaving  the  laser,  the  beam  was  directed  through  a  half-wave  plate  that 
served  as  a  beam  attenuator  when  necessary.  Following  the  half-wave  plate  was  an 
apparatus  for  isolation  of  the  laser  from  the  phase-conjugate  beam.  The  apparatus 
used  was  a  faraday  rotator,  manufactured  by  Optics  For  Research  Inc.,  consisting  of  a 
25mm  long  rod  of  Hoya  FR5  glass  surrounded  by  a  set  of  permanent  magnets.  The 
isolator  was  set  for  optimal  rotation  of  the  plane  of  polarization  of  the  beam  at  45° 
and  placed  between  twn  polarizers  with  axes  oriented  at  45°  to  each  other.  This 
configuration  gave  a  backward  attenuation  of  approximately  500.  Without  the  isolator 
the  phase-conjugate  beam  entered  the  laser  cavity  and  caused  large  fluctuations  in 
power  output  as  well  as  mode  hopping.  The  polarization  of  the  beam  upon  leaving  the 
isolator  was  such  that  it  would  enter  the  crystal  extraordinarily  polarized,  which  is 
required  for  the  formation  of  a  SPPCM. 

After  passing  through  the  isolator  the  beam  was  expanded  by  a  factor  of  6.25 
with  an  expanding  telescope.  This  allowed  the  insenion  of  transparencies  into  the 
beam  when  desired,  enabling  the  verification  that  an  observed  beam  is  a  phase- 
conjugate  beam.  .Also  included  in  the  telescope  apparatus  was  a  0.05  mm  pinhole  that 
could  be  removed  at  will.  This  pinhole  provided  the  ability  to  make  observations  with 
a  Gaussian  beam  (pinhole  inserted)  or  with  what  was  eventually  termed  a  mottled 
Gaussian  beam  (wthout  the  pinhole).  The  mottled  Gaussian  was  so  named  because, 
although  the  beam  without  the  pinhole  was  generally  Gaussian,  it  had  a  mottled 

no  '.'.J  'n.e  'c.am  '.v.is  iircctou  ’.nrough  a  beamsplitter  xnth  'o"'-'  • 

transnussion  at  514  nm.  This  beamsplitter  not  only  provided  the  means  for  viewing  the 
phase-comugate  signal  as  it  returned  along  the  incident  path,  but  also  provided  a 
reference  beam  that  was  used  to  look  for  frequency  shifts  of  the  phase-conjugate  beam. 
To  .ock  .‘'or  a  frequency  shift  the  phase-conjugate  return  was  reflected  from  the 
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Figure  A.  1  Experimental  Apparatus. 


b€amsplutcr  and  redirected  to  interfere  with  a  portion  of  the  inadent  beam  reflected 
from  the  front  side  of  the  beamsphtier.  The  resultmg  interference  pattern  was  then 
observed  for  evidence  of  movement. 

Upon  passage  through  the  beamsplitter,  the  beam,  with  power  of  75  ±5  mW, 
was  erher  mcdif.ed  by  a  lens,  slit,  k.nife-edge  or  a  combination  thereof,  or  simply 
collimated  and  directed  toward  the  crystal.  As  the  phase-conjugate  beam  returned 
along  the  incident  path,  the  beamsplitter  reflected  a  portion  of  it  and  directed  it  to  a 
large  area  silicon  PIN  photodiode  reverse  biased  at  22  volts  and  terminated  with  a  100 
ohm  resistor.  The  output  of  the  detector  was  sent  to  a  digital  millivolt  meter,  an 
oscilloscope  and  a  chart  recorder.  When  an  image  was  placed  in  the  beam  the  detector 
was  removed  and  replaced  wnth  a  white  screen  for  viewing.  When  frequency  shift 
information  was  desired,  the  detector  was  replaced  with  a  mirror  and  the  two  beams 
were  directed  toward  a  white  screen  for  observation  of  the  resulting  interference 
pattern 

The  cry  stal  was  any  one  of  four  nearly  cubic  single  cry  stal  of  BaTiO^  measuring 
5±'f2  mm  on  a  side.  The  sides  of  the  crystals  were  polished,  but  it  is  im.porant  to 
note  that  the  sides  were  not  exactly  the  same  size  and  therefore  the  sides  were  not 
exactly  orthogonal.  The  crystals  were  mounted  on  a  platform  that  could  be  raised, 
lowered,  rotated  and  tilted.  Observations  of  internal  beams  were  made  visually  from 
above  and  were  recorded  by  film  or  video  tape  Both  the  still  and  video  cameras  used  a 
f  ■  3  5mm  Cannon  FL  lens  mounted  on  a  single  track  bellows  extended  to  provide  a 
full  fra.me  view  of  the  crystal. 

Mans  of  the  photographs  presented  here  were  taken  from  the  video  tape  directly 
oT  c*  1  ’cie’.  >i>'n  screen  using  a  bx~  cm  format  still  camera  The  desired  frame  of  the 
video  was  digi'i/ed  and  placed  on  a  hard  disk  .''or  recall  and  photographing,  the 
numbers  dispiavcd  in  'he  photographs  are  a  result  of  the  equipment  used  to  do  this 
‘‘inu.  pbo’ograrh'.  ind  ire  nor  pertinent  to  the  investigation  When  this  method  was 
...Ned,  wis  pussiOie  'o  m.jKC  negative  images  of  the  desired  photographs,  which 
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±  3® 0.  In  the  photographs  presented,  except  where  noted  otherwise,  the  incident  beam 
enters  vertically  from  the  top  and  is  visible  only  upon  entering  the  crystal. 

3.  RESULTS  AND  ANALYSIS  OF  DATA 
a.  Experimental  variables 

The  number  of  variables  involved  in  SPPC  is  immense,  however  1  have  tned  to 
concentrate  on  those  which  have  emerged  as  the  most  important.  The  most  important 
factors  in  the  production  of  various  beams  differ,  but  they  generally  may  be  considered 
as  dependent  on  beam  profile  or  the  geometry  of  the  beam-crystal  interface. 

Lenses,  knife-edges,  slits  and  pinholes  were  used  to  change  the  profile  of  the 
incident  beam  prior  to  the  beam  entenng  the  ciy  stal.  The  effect  of  these  modifications 
of  the  incident  beam  is  especially  obvious  when  the  beam  profile  has  a  large  gradient, 
such  as  when  a  knife-edge  is  inserted  into  the  beam,  since  this  creates  a  very  clear 
delineation  in  the  lighted  and  unlighted  areas  of  the  cry  stal  and  introduces  new  wave 
vectors  into  the  beam.  The  effect  of  the  lens  to  crystal  distance  on  phase-conjugate 
reflectivity  has  been  studied  by  A.  V.  Nowak  and  is  still  not  well  understood  [Ref  37], 

The  geometry  of  the  beam-crystal  Interface  is  probably  the  most  important  of 
all  variables  in  the  system.  Changes  in  angle  of  less  than  one  degree  and  lateral  shifts 
in  the  beam  of  less  than  0.5  mm  can  cause  drastic  differences  in  internal  beams 
configurations.  The  following  examples  are  presented  to  indicate  the  drastic  effects 
observed  with  these  type  of  changes. 

The  configuration  shown  in  Figure  A. 2  is  a  typical  experimental  configuration 
that  produced  a  phase-conjugate  signal  from  the  BaTiOj  ciystal.  The  phase  conjugate 
reflectivity  was  33'’  o.  Note  the  fanned  beams  at  the  lower  left  comer  of  the  ciy  stal. 

Figure  .A. 3  demonstrates  the  extreme  effects  made  possible  when  changing  the 
beam  profile  by  placing  a  knife-edge  m  the  incident  beam.  Unless  otherwise  indicated, 
the  knife-edge  always  intercepted  the  beam  on  the  right  side  and  was  placed  behind  any 
lenses  Note  the  extreme  curvature  of  the  internal  beam  observed  when  blocking  hall' 
the  beam  \sTth  a  knife-edge  This  also  reduced  the  phase-conjugate  refiection  to  7ero 

\  ,argc  eTect  is  also  seen  when  relocating  the  beam  entrance  position  b\  a 
iaterii  siuf:  m  the  bcj.m-cr. sta!  intersection  As  the  beam  is  moved  across  the  l.ise  of 
'ne  cr.sta..  the  ph  ise-ccnjuga'e  rel'.ectivuy  will  change,  beginning  with  no  phase 
ccn;uga’i;n  it  one  end  and  ending  with  the  same  result  at  the  other  In  between  the 


Figure  A  2  Normal  configuration  for  BaTiOj  exhibiting  SPPC. 


two  extremes,  the  phase-conjugate  reflectivity  may  range  as  high  as  6<)'’o  or  be 
nonexistent  depending  on  other  factors.  Similarly,  a  change  in  the  angle  of  the  beam 
entering  the  crystal  wili  have  a  drastic  effect  on  the  phase-conjugate  reflectiuty  and 
beam  pattern.  It  is  generally  unpredictable  and  small  changes  in  angle  (  ha'-e 
^een  observed  to  have  no  effect,  to  destroy  a  pattern  and  eliminate  all  phase-conjugate 
reflectivity,  or  to  increase  the  phase-conjugate  reflectivity  as  much  as  SCO”  ), 
b  Comple.x  Au.\iliary  Beam  Patterns 

/.  The  Total  Internal  Reflection  {TIR)  Ring 

One  of  the  most  interesting  effects  discovered  was  the  TIR  ring  Figure 
.A  -1  shows  the  first  example,  to  my  knowledge,  of  a  totally  internally  pumped  ring 
IS. on  'Vi'hm  i  BiTiO.  .r  <’ai-  e:Te;.t  ''•a'  mh.e'.c  ;  ..s,;!,:  i 

,  :  ...  !  .  Vc  ;  j  i  '  .j,  .  ■■c  .  . ^ 

■cr.-.,  :  -.m  irorn  tne  lace  oi  tne  crj.stai.  1  he  beam  was  incident  at  -id"  to  tb.e  axis  o; 

the  crystal  and  the  knife-edge  was  placed  on  the  side  of  the  beam  nearest  tine 
The  effect  was  occasionally  obtained  -.vhiie  obscuring  the  other  hall  of  the  beam,  and 
'.vas  once  obtained  w-.th  a  ver’icailv  ahcn.ed  slit  A  V  Nowak  once  obtained  tlie  rme 


Figure  A. 3  EfTect  of  placing  a  knife  edge  in  the  incident  beam. 

with  ;ust  the  focusing  lens,  indicating  that  the  TIR  ring  may  be  dependent  more  on  the 
availabiiits  of  light  scattered  at  a  certain  angle  than  on  the  incident  beam  profile. 

In  the  case  shown  in  Figure  A. 4  the  phase-conjugate  signal  took  about  '<• 
seconds  to  appear.  The  beam  within  the  cry'stal  slowly  moved  to  the  corner  until  the 
configuration  of  the  beam  arrived  as  shown,  but  without  the  wng.  Then.  siowK,  the 
ring  appeared  o\cr  a  three  second  interval.  the  ring  appeared,  the  phase  conjugare 
reilecti'. itv  dropped  from  to  b"'’.  ,A  drop  in  phase-conjugate  refectr.  .'v  Aa> 

observeu  ;n  all  examples  of  ring  formation  for  those  cases  where  the  lefecti'-itv  uas 
measured  In  some  cases,  several  nearly  overlapping  rings  were  observed,  as  m  1  icurc 
.\  4  The  ring  in  Figure  .-X  4  -.vas  steady,  fading  cccasionalh  ion  the  order  o;  lue 
e'.er',  ■  .e  m^nu’es.  '•'c*  hv.  i’.  'e' ..r’lin^  '■  v  "'ct  '\er:'  ■ 

iiiii  :rx  :j;..t.on  was  unuer  i  tnen  tnc  ring  '.vouiu  reappear  lians.u’ioii  a, '-.c 

tne  v.r'.snh  on  ':.e  order  c!  n  '  mm  in  either  direct. on  did  not  disturb  the  nnc. 

Tans. anon  .vou.l  destrov  it  I  or  other  ccnf gurations,  the  ring  aas  \  pieser.eu  !  ■ 
one  d.rect.on  o!  cr. stal  mot. .or.  ■  «ee  Secti  m  4i 
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Figure  A. 4  The  TIR  Ring. 

The  direction  of  circulation  within  the  rings  was  determined  b\  extracting 
some  of  the  light  out  of  the  ring  into  a  second  crxstal  of  BaTiOj  After  producing  the 
ring  in  one  crcstai.  a  drop  of  index  matching  liquid  imethvlene  iodide  CM,!,’  diluted 
1  1  uith  2-propanoi.  n  "  1  6)  was  placed  onto  the  surface  of  a  second  cnstal  and  ’lie 
second  cTV'tal  .cas  moved  into  contact  with  the  first  Since  the  loss  into  the  second 
.r'>st.*.  Was  too  great  to  aiiovv  the  continuation  of  the  ring,  the  operation  was  rcc  ided 
'he-  . :  Jeo  .a.'nera  and  the  frame  of  the  .nstant  'aas  ot  contact  was  stud.cu  a'  i  .a-  rr 
'.me  Stjd-.  of  the  video  revealed  that  the  light  within  the  ring  could  pr^  paga'e 
.o',.r.te:'c.ocKWisc.  or  both  wavs  simuitaneouslv  depending  on  the  parr.c-.  ar  •  -ig 
Hi  '.vc'.  cr,  .n  'he  majoritc  of  cases,  oniv  'he  vOuntcrciocK'.v ise  direction  aav  '''"crvc. 
f  c  \  '  -lea''',  'h  's  'h.s  .  ,,.n'er  '  sr  i-  ,  v  x.;-’  • 


1  ne  ring  aittnn  the  .r.stai  .s  pr..)hao.\  produced  r^\  twon'cirn  ^ 
’^e’ween  'he  -nain  hcarn  and  scat'ered  .ight  that  linJs  a  dosed  l  ''>p  ..a 
'Ctlectuns  aroji'.J  'he  crvstal  Ih.s  somew.hat  i  ss\  IIR  ring  takes  pauc  ’• 
.n^ident  ^eam,  -.chich  expia.ns  :he  ,ir''p  in  p'lasc-d. npugate  on'lcc'.v  .'v  'oe.-', 
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Figure  A  6  The  Diamond  configuration 

<harp  diamond  appears  in  tlic  crystal  the  phase-conjugate  signal  increased  considcrabU 
I  to  Ml  tirr.CN  •'■hat  it  is  just  prior  to  its  appearance)  The  sharp  diamond  probihh 
l^  .1  specul.ir  rctlection  if  the  phase-(.onjugate  signal,  uhich  reflects  back  to  ttic 
.:x:dent  ‘'earn  internal  reflection  Sin«.c  a  second  point  for  \1  is  created,  one 
a  j  .Id  expcc'  ti.c  obser'.ed  increase  in  the  pfiasc-conjugate  signal  upon  the  appea.'ancc 
’.'.c  diamond 

A  :asv.,r.at.ng  and  unexpia.ned  ■.  anat  .or.  o;  the  sharp  diamond  is  the  dill  use 
d.irr.ond  The  diifjse  diamond,  sho'-sn  m  Figure  .A  begins  .n  a  Lonf'igurat.mt  similar 
•  ■  ".e  d.ammd  shi'-.sn  ;n  Figure  A  ^ut  l^  tfie  ^rss-ai  is  lateralis  rnosed,  a  portion  of 
me  frt'.t  '''earn  .  ;mpO'ing  the  sharp  diamond  s'ass  to  .rcate  a  diamond  :n  the  nchr 
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Figure  A. 7  The  Diffuse  Diamond  configuration. 

reflection  of  the  input  beam  within  the  crystal  creating  a  diamond  but  there  is  no 
reflection  of  the  sccondar>'  beam  off  of  the  rear  face  of  the  crystal,  Absorption  of  the 
reficcted  beam  by  two-beam  coupling  with  the  incident  beam  within  the  first  few 
nucrometers  of  the  crvstal  edge  is  probably  responsible  for  the  lack  of  specular 
re.fcction  in  this  and  other  similar  cases.  There  is  a  similar  theory  which  account's  for 
'he  lack  of  specular  rcficction  of  the  phase-conjugate  b^"  im  off  of  surfaces  such  as  the 
input  face  of  tne  cry  stal.  Pepper  (Ref.  56|  proposes  that  the  lack  of  a  retlection  of  the 
phase-^cnjugate  beam  is  due  to  destructive  interference  between  the  incident  beam 
refection,  and  a  beam  created  by  the  phase  conjugation  of  the  retlection  of  the  original 
phase- :  mjugate  beam.  This  has  been  borne  out  in  c.xperimcnt.  However,  as  F  igure 
\  -i,‘m.::ns'ra'es.  i’  cer’ain  ancies  'he  '•'e^'.f.m  'ct'.ec'ion  ;s  su.Ticientb.  ■  ' 

. .  _ _  .i.cr'ta.  i;  ■■  .:i .  'na.'.  i"-  /:  c- 

iac\  o:  a  retlection  iiave  yet  to  be  fuiiy  determined 
4.  Diagonals 

r.gurc  .-X  111  illustrates  the  diagonals  The  diagonals  faded  in  and  cut  i'-  '..e 
ph.ise-LonjUga'c  reilectr-itv  oscillated  from  12’.,  to  21’"  with  a  period  on  the  crder  n 


Figure  A. 8  Example  of  lack  of  reflection  of  the  incident  beam 

%vnthin  the  cr>$tal. 


T;  e  reflectivjrv  was  strongest  when  the  diagonals  were  brightest  When  the 
incidence  was  increased  to  06*^  a  similar  configuration  appeared,  but  wuh  the 
jnjug.ite  refectivtty  oscillating  from  one  to  .and  specuiar  refection  olf  of 
s  ;.i^e  lading  as  the  diagonals  became  brighter.  .At  another  angle  of  inc.dcnco 
. or.l.guraticn  (  62^’’  it  was  noted  that  increasing  the  power  of  the  inpin  beam 
nun  ;:e  a.l  oscillations  ot  the  pha'ic-ccnjugate  refection  while  elinunann^  the 
If,;>  -or.i'iguration  hac  also  been  seen  vs'ithout  anv  mcasura^'ie  r'ha^c- 
*  h'  r  0 ! c  w ! r* 

Temporal  Effects  and  Frequencs  Shifts 
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Figure  A. 9  Example  of  multiple  reflections  within  the  cr^’stal. 

/.  Amplitude  Oscillations 

Frequently,  amplitude  oscillations  were  observed  in  the  phase-conjugate 
signal  of  many  of  the  previously  discussed  examples.  These  include  Figures  .A. 2,  .A  d, 
.A.b,  .A.S.  .A. 9.  The  oscillations,  probably  similar  to  those  reported  by  Valley  and 
Dunning  in  an  external  resonating  cavity  iRef  ?9].  were  most  obvious  wlion  a  mottled 
':^eam  was  focused  into  the  sample,  but  they  also  appeared  in  some  conf.guraticns  wnh 
a.i  "earn  profiles.  The  observed  oscillations  often  appeared  to  he  chaotic,  and  mu^f. 
'.vork,  lias  been  reported  recently  concerning  the  observation  of  chaotic  oscillations  ,:i 
BaT.O,  Ref  5d,35]  However,  some  of  the  observed  oscillations  were  \erx  near!', 
per; 'die  ■  per.odic  •.'.••.thin  I'f’  i  of  some  mean  period'  with  penods  ranging  from  'o 


a;’:.  1  fwuscd  Ga'v^ssian  beam  incident  'n  the  center  of  the  crystal  do''  to  ’ne 
I'ljNC  .  vw.datmns  appeared  ;mm.edi...telv  -upon  the  initiation  of  phase  ..■I'n  ugat  :. 


Ihesc  •.'•ere  m.untaineu  as  sho'»sn  lor  a  period  .;l  anout  one  nunutc  and  tlien  i‘'r',.p' 
i-hanc-od  t  .'  ti.e  osci...iti..'ns  sha-.cr.  Fieure  \  12  'A'lthin  a  t'e'.s  tninutcs  •lie  rie  ;; 


Figure  A.  10  The  Diagonals  configuration 

oscillations  reappeared  Fast  oscillations  ( 10  Hz*  amplitude  modulated  bs  sloucr 
oscillav.oas  '  "05  Hz',  and  oscillations  with  faster  beats  superimposed  on  them  '.\e:e 
u.so  observed  Nowak,  has  shown  tnat  the  frequenev  of  these  amplitude  oscillations  ;s 
proportional  to  the  intensits  of  the  incident  beam  jRef  o';  and  this  has  been  ccntlrmed 
■'V  ot.ners  Ref  3"' 

These  amplitude  osculations  of  the  phasc-con;ugatc  signal  are  probabh.  a 
result  c:  frequcncs  shifts  of  those  auxilian  '■'earns  that  acvount  fo-  the  fo-mation  .'i  ;:ie 
grat.r.c  Frequenev  sh.fts  of  the  auxiharv  beams  vsere  proposed  b\  I  emberg  and 
nuw.ner  Ref  J''  to  explain  observed  irequcr.cx  sweeping  in  the  output  o:  a  ring 
r'.’s  •:  It  or  using  B-iIit),  and  in  in  .niso.ated  A:'  laser  producing  phase  vOn  ..g.r,  on 
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Figure  A.  1 2  Amplitude  oscillations  which  progressed  from 
those  of  Figure  A.  1 1. 

:c'*e;at;on  of  amplitude  and  frequency  shift  of  the  phase-conjugate  beam  may  have 
^eer  exp! lined  bx  MacDonald  and  Feinberg  when  they  showed  that  the  four-wave 
muxing  signal  falls  off  smoothly  as  the  frequency  difference  between  the  reflected  and 
;rc;dent  beam,  becomes  larger  than  ~  1  Hz.  Although  the  mechanism  for  this 
:re.;..cncv  sh..''ting  is  still  under  investigation,  Lam  [Ref  3]  has  proposed  that  the  bulk 
r.nc:x voltaic  effect  gives  rise  to  an  intensity-dependent  frequency  shift.  This  implies 
•nat  tne  frequency  is  dependent  on  the  amplitude  instead  of  vice  versa.  .Another 
exn'anaticn  may  be  that  there  is  consistent  movement  of  the  grating  when  the  reflected 
r-n.p  vave  form.s  a  loop  (see  Figure  .A. 2  for  example).  In  this  case  each  circuit  around 
•fs  l  .'cp  -.v-.;:  Dcpnler-chift  the  pumpwave  and  the  observed  effect  may  be  expected. 

'  ■.ere  'r.."  ;ne  'n^crvou  nstar.ce  vnere  :ne  frequencv  i;;feren..e 

e  ee;.  ...e  . .. c  oeair.')  .'as  non-zero  ana  constant  i^no  photograph  available).  The 
pn  isc-'. 'nyugatc  return  was  frequency-shifted  from  the  input  beam,  but  there  was  no 
v.^eeping  T'r.at  is.  there  was  obv-.ous  movement  of  the  interference  fringes  at  a 
r.siant  ve.ccitv  This  effect  resulted  from  a  mottled  Gaussian  beam  incident  on  the 


I 


cn’stal  at  an  angle  of  35°  to  the  c-axis.  Hcvvever  this  .s  ar.  jronia.  i 

yet  to  be  repeated  or  explained. 

3.  Influence  of  Returning  the  transmitted  Beam 

This  effect  is  similar  to  that  reported  by  f  embcrg  !Rc:  2"  i  .  v  ; 
discovered  while  trying  to  stabilize  the  amplitude  oscillations  :!  the  r- 
signal  discussed  above.  The  effect  was  seen  using  an  arrangement  ’  - 

nonperiodic  amplitude  oscillations  of  the  phase-conjugate  rcHccticn  rangir  c  -  ■  .' 

to  46%.  In  this  configuration  a  100%  reflectance  mirror  was  placed  behind  -r  e  : 
directing  the  transmitted  beam  back  through  the  crystal,  but  net  direi.-.'  ,  ■ 

incident  path.  All  oscillations  stopped,  leaving  the  phase-conjugate  ref'.e.n  n 

10°'o. 

The  explanation  for  this  effect  is  probably  found  in  the  i-ompc-  •.  • 
between  beam  patterns.  In  any  instance,  the  observed  beam  pattern  is  the  or  e 
exhibits  the  maximum  gain.  When  oscillation  occurs  between  two  compet;:  g  ''e,.:-. 
patterns  it  is  apparent  that,  for  some  reason,  the  gain  of  one  is  falling  abuse  and  -r.er. 
below  the  possible  gain  of  another.  The  input  of  a  beam  of  light  from  any  sou.-^ve  '.s.';. 
brightness  on  the  order  of  the  internal  beams  will  cause  a  change  in  the  gam  .'!  ',,e 
system  due  to  some  or  all  of  the  possible  nonlinear  effects.  In  the  case  of  returning  th.e 
transmitted  beam,  it  appears  that  the  gain  of  one  of  the  competing  patterns  is 
sufficiently  reduced  to  eliminate  it  from  the  possible  patterns, 
d.  The  Free  Standing  Auxiliary  Beam  Pattern 

This  section  presents  an  occasionally  observed  result  due  to  a  small  crystal 
motion.  Here  the  crystal  was  translated  at  right  angles  to  the  beam  path  while  the 
internal  beam  pattern  was  carefully  observed  for  any  changes. 

Note  the  beam  configuration  seen  in  Figure  A.  13.  When  the  crystal  was 
shifted  slightly  to  the  left,  the  auxiliary  beam  pattern  moved  with  the  crystal  as  shown 
in  Figure  A.  14  instead  of  changing  to  adapt  to  the  new  crystal  orientation.  Note, 
however,  that  some  slight  distortion  of  the  pattern  did  occur.  In  this  case  the  phase- 
conjugate  re.'lectivity  was  originally  32%  '«vith  small  amplitude  oscillations  on  the  order 
of  5''d.  There  'vas  no  evidence  of  Treuuencv  -.niiting  of  the  ohase-conida  i.-‘  'c 
Upon  translation  of  the  crystal,  the  phase-conjugate  reflectivity  diminished  to  23’  > 
with  oscillations  of  the  same  amplitude  as  before  the  translation. 

After  translation,  any  rotational  motion  of  the  crystal  destroys  the  auxiliary 
beam  and  a  new  one  appears  that  is  similar  to  the  original  pattern  (i.e.  Figure  .X  L') 


Figure  A  14  Free  standing  beams  created  by  translating  the 
crystal  '^0.5min  from  the  configuration  of  Figure  A.  13. 
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10 

20 

30 

40 

50 

60 

70 

80 

90 

100 

110 

120 

130 

140 

150 

160 

170 

180 

190 

200 

210 

220 

230 

240 

250 

260 

270 

280 

290 

300 

310 

320 

330 

340 

350 

360 

370 

380 

390 

400 

410 

420 

430 

440 

450 

460 

430 

430 

500 

510 

520 

530 

540 

550 

560 

570 

580 


1PSIM_3 

I  A  SIMULATION  OF  INTERNAL  BEAM  RESPONSE  IN  BARIUM  TITANATE 

!  IBEAM  3  SIZE  OF  INPUT  BEAM 

I  IBEAMF  »  ACTUAL  NUMBER  OF  POINTS  IN  WAVEFRONT  AFTER 
I  DESIGNATION  OF  ANGLE  PHI 

I  ISIZE  «  SIZE  OF  CRYSTAL 
!  ICRYS  0  »  CRYSTAL  GRID 

!  KEEPO  »  ARRAY  TO  KEEP  TRACK  OF  CURRENTLY  OCCUPIED  SPACES  IN  THE 
!  GRID,  2=COL,  l=ROW 

!  IPHl  ■  INPUT  ANGLE  OF  INCIDENCE 
I  PHI  ■  ANGLE  OF  INCIDENT  BEAM  (REAL) 

I  PHIC  ■  ANGLE  OF  INCIDENT  BEAM  AFTER  ENTERING  CRYSTAL 
!  IX  a  NUMBER  OF  COLUMN  POINTS  ACROSS  WAVEFRONT 
lY  =  NUMBER  OF  ROW  POINTS  ACROSS  WAVEFRONT 
!  ADD()  a  AMOUNT  TO  ADD  TO  Y  FOR  EACH  X  MOVED 
!  I ROW  a  ROW  NUMBER 

I  I COL  a  COLUMN  NUMBER 

!  ITOP  a  highest  index  OF  REFRACTION 
I  IBOT  a  LOWEST  INDEX  OF  REFRACTION 

!  GRAD  a  GRADIENT  OF  INDEX  OF  REFRACTION  ACROSS  THE  BEAM 
!  CRYS  a  real  value  OF  INDEX  OF  REFRACTION  (BEFORE  FIXING) 

!  INDEX  a  INDEX  OF  REFRACTION  OF  LOCATION 

REALX  a  PLACE  FOR  KEEPING  TRACK  OF  AMOUNT  TO  ADD  TO  COLUMN 
(OR  ROW  FOR  REALY)  WHEN  MOVING  DOWN  THE  BEAM 
THIS  ENSURES  THAT  ALL  FRACTIONS  OF  A  GRID  SPACE 
ARE  NOT  LOST 
REALY  a  SEE  REALX 

REALAD  0  a  SAME  AS  REALX  ABOVE  BUT  ONE  SPACE  FOR  EACH  POINT. 

USED  DURING  ACTUAL  WAVEFRONT  MOVEMENT 
IPERMCO  a  STARTING  COLUMN  FOR  WAVEFRONT  POINT 

IMARKB  aMARKER  TO  KNOW  WHEN  REFLECTION  OFF  OF  BOTTOM  HAS  OCCURRED 
IMARKS  a  marker  TO  KNOW  WHEN  REFLECTION  OFF  OF  SIDE  HAS  OCCURRED 
ICHECK  a  index  of  TIME  PERIODS 
RSTOP  a  number  of  TIME  PERIODS  IN  SIMULATION 


IDIR  a  number  of  time  PERIODS  BETWEEN  COMPUTATIONS 
OF  DIRECTION  OF  WAVEFRONT 
MDIR  a  COUNTER  FOR  CHECKING  AGAINST  IDIR 
XX  a  HORIZONTAL  DISTANCE  BETWEEN  END  POINTS  OF  WAVEFRONT 

w  = 


;  :hik 


VERTICAL  DISTANCE  BETWEEN  END  POINTS  OF  WAVEFRONT 
=  VARIABLE  VSED  TC  IHANGE  DIRECTION  AFTER  REFLECTION 
.IPJiC;,  =  ''.AIFI/:  -JSED  TO  :L\RK  the  incident  beam  path  30  THAT 
THE  INDICES  SET  UP  3Y  THE  INCIDENT  BEAM  DOMINATE 
MEAN  a  MEAN  VALUE  OF  THE  INDEX  OF  REFRACTION 
REF  a  actual  index  OF  REFRACTION  OF  THE  PHOTOREFRACTOR  BEING 
SIMULATED 

IPLACE  a  POINT  THAT  BEAM  STARTS  IN  CRYSTAL 

IGRAPH  a  marker  FOR  POINTS  TO  BE  PLOTTED  (SET  FOR  100  POINTS) 
IGRMAX  a  NUMBER  OF  ITERATIONS  BETWEEN  PLOTTING  WAVEFRONTS 
CAXIS  a  FRACTION  OF  TOTAL  INDEX  CHANGE  DUE  TO  ANGLE  BEAM 
MAKES  WITH  THEE  C-AXIS  IN  THE  CRYSTAL 
CINDEX  a  REAL  VALUE  OF  INDEX  OF  REFRACTION  USED  TO  COMPUTE 
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590 

600 

610 

620 

630 

640 

640 

650 

660 

670 

680 

690 

700 

710 

720 

730 

740 

750 

760 

770 

780 

790 

800 

810 

820 

830 

840 

850 

850 

870 

880 

890 

900 

910 

920 

930 

940 

950 

960 

970 

980 

990 

1000 

1010 

1020 

1030 

1040 

1050 

1060 

1070 

1080 

1090 

1100 

1120 

1130 

1140 

1150 

.150 

.1"0 

>13C 

1190 

1200 

1210 

1220 

1230 

1240 

1250 

1260 

1270 

1280 


I  THE  CHANGE  DUE  TO  THE  ANGLE  TO  C-AXIS 

!  lOUTO  ■  MATRIX  FOR  HOLDING  VALUES  OF  PLOTTED  POINTS 

1  ICOUT  »  COUNTER  FOR  PLACE  OF  POINT  IN  OUTPUT  FILE 

I  ISKIP  »  MARKER  TO  CHECK  IF  WAVEFRONT  IS  UNDERGOING  REFLECTION 
I  REFLCHKB  -MARKER  TO  KNOW  WHEN  REFLECTION  OF  WAVEFRONT  STARTS/ENDS 
I  REFLCHKS  »  SAME  BUT  FOR  SIDE 

I  OPl  »  NUMBER  OF  POINTS  IN  WAVEFRONT  TO  LOOK  AT  FOR  RECALCULATION 
I  THIS  EQUALS  IBEAMF/PIECE 

!  ISETMARK  =  MARKER  TO  TELL  WHEN  TO  STOP  LOOKING  FOR  REFLECTION 

!  PIECE  =  NUMBER  OF  PIECES  OF  WAVEFRONT 

!  ADDYO  *  FRACTION  OF  ONE  BOX  THAT  EACH  POINT  STARTS  OUT  AHEAD  OF 

!  ITS  ROW 

1 
! 
i 


INTEGER  Ibeam,  Ibeamf ,Isize,Icrys(500,500) 

INTEGER  Ke«  ( 2 . 50 ) ,  Iphi ,  Ix ,  ly ,  I  row ,  Icol ,  I  top ,  Ibo  t ,  Re  f  Ichks  ( 50 ) 
INTEGER  Reflchkb(50) ,  Ikfac, This, That 

INTEGER  Imdex, Imaexp(50) , Ipermc(50) , Imarkb(50) , Itime(50) , Iskipb 
INTEGER  Imarks ( 50 ) , lair ,Mdir , Ichng, Icout , Iskips , I temp , Opl , Piece 
INTEGER  Itrac( 500 , 500) .Mean, Iplace , Igraph, loutTl , 50,101) , Isetmark 
REAL  Phi,Phic,Add(50) ,Grad,Crys,Realx,Realad(50) ,Rstop,Grmax 
REAL  Xx , Yy , Re t , Caxis ( 50 ) , Cindex , Rcheck , Kf ac , Realy , Addy { 50 ) 


Iskipb 

.Piece 


! 


!  NOTE  THAT  IF  BEAM  SIZE  EXCEEDS  50  YOU  MUST  REDIMENSION 
! 

ON  ERROR  GOSUB  5430  t  WRITES  FILE  TO  DISK  IF  ERROR 
PRINTER  IS  1  !  SET  TERMINAL  OUTPUT 

GRAPHICS  OFF 

OUTPUT  KBD;  K;  'CLEARS  SCREEN 

ON  KEY  2  LABEL  SAVE  &STOP  GOSUB  5430 
PRINT  TO  STOP  PROGRAM  AND  SAVE  DATA  PRESS  f2 


INPUT  FACTORS 


INPUT  NAME  OF  DATA  FILE  TO  BE  USED  ,Name$ 

PRINT  DATA  FILE:  ,Name$ 

INPUT  SIZE  OF  CRYSTAL  TO  BE  USED  ,Isize 
PRINT  CRYSTAL  SIZE  IS  ,Isize,  ON  A  SIDE 
INPUT  BEAM  SIZE  IS  , Ibeam 
PRINT  WITH  A  BEAM  OF  SIZE  ,  Ibeam 

INPUT  PLACE  ON  CRYSTAL  FOR  BEAM  TO  ENTER  IS  ,  Ip lace 

PRINT  ENTERING  AT  ,  Iplace 

INPUT  HIGH  VALUE  OF  INDEX  IS  ,  I top 

PRINT  LOW  VALUE  OF  INDEX  IS  ,Ibot 

PRINT  INDICES  RANGE  FROM  ,Itop,  TO  , Ibot 

INPUT  INITIAL  ANGLE  IN  DEGREES  IS  (89  DEG  MAX)  ,  Iphi 

T1=TIMEDATE 

I 

I 

:NI7:.\L;ZE  -.ikRIABLES 

.^checx  =  0 
Mdir=0 
Ref =2. 4 
lcout=0 
lsetmark=0 

Piece=5  !  REM  TO  CHANGE  KK  IF  YOU  CHANGE  PIECE 


86 


1290 

1300 

1310 

1320 

1330 

1340 

1350 

1360 

1370 

1380 

1390 

1400 

1410 

1420 

1430 

1440 

1450 

1460 

1470 

1480 

1490 

1500 

1510 

1520 

1530 

1540 

1550 

1560 

1570 

1580 

1590 

1600 

1610 

1620 

1630 

1640 

1650 

1660 

1670 

1680 

1690 

1700 

1710 

1720 

1730 

1740 

1750 

1760 

1770 

1780 

1790 

1800 

1810 

1820 

1830 

1840 

1350 

1360 

.37C 

1890 

1900 

1910 

1920 

1930 

1940 

1950 

1960 

1970 

1980 


PRINT  NUMBER  OF  PIECES  IN  WAVEFRONT  ,  Piece 

I 

! 


f 

i 

I 


PUT  FACTORS  IN  OUTPUT  FILE 


lout(0,0,0)«lsize 
lout(0,l,0)«lbeaffl 
loutil ,0,0)»lplace 
Iout{ 1 , 1 ,0 )«Itop 
Iout{0,0.l6l)»Ibot 
lout(0, 1, 101}*Iphi 
lout ( 1 , 0 , 101 )«Piece 
lgr«ph«0 

Ratop-Itop*2  !SET  TIME  TOTAL  PERIODS 

Rstop»Rstop*Isize  I  USE  TWO  LINES  TO  KEEP  FROM  INTEGER  PROBLEMS 
PRINT  TOTAL  TIME  »  ,  Rstop 

Grmex-Rstop/lOO  I  SET  FOR  100  WAVEFRONTS  PLOTTED 
GOSUB  5170  ISET  UP  GRAPHICS  DISPLAY 
I 

I  TIME  BETWEEN  NONLINEAR  CALCULATIONS 
Idir*5*Itop 
! 

! 

I  INITIALIZE  ARRAY 
! 

FOR  1*1  TO  Ibeam 


>(1,I)»0 


Keep 
Keep(  . 

Itime(I)=0 
Refichks(l)“0 

ReflchkbU)=0 

NEXT  I 
! 

! 

! 

!  CONVERT  PHI  FROM  DEGREES  TO  RADIANS  AND  FIND  ANGLE  INSIDE  CRYSTAL 
Phi»Iphi 

Phi*Phi*(2*PI/360) 

Phic*(SIN(Phi)/Ref) 

Phic*ASN(Phic) 


CHANGE  SIZE  OF  BEAM  DEPENDENT  ON  THE  INPUT  ANGLE 


X=Ibeam*COS (Phic ) 

Y=Ibeam*SIN(Phic) 

Iy=Y 

Ix=X 

Ibeamf=Ix 

Opl=Ibeamf/Piece  IISIZE  OF  PIECE  TO  LOOK  AT  WHEN  CALCULATING  DIRECTION 

Add(l)=Y/X 

Addl*Add(l) 

IF  Iplace>(Isize-Ibeanif)  THEN  I  MAKE  SURE  IPLACE  IS  NOT  TOO  LARGE 
Iplace  *Isize-Ibeamf 
END  IF 


! 

! 

! 

Irow*0 

Realy=0 

Icol*Iplace 


FILE  WAVEFRONT  COORD 
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1990 

2000 

2010 

2020 

2030 

2040 

2050 

2060 

2070 

2080 

2090 

2100 

2110 

2120 

2130 

2140 

2150 

2160 

2170 

2180 

2190 

2200 

2210 

2220 

2230 

2240 

2250 

2260 

2270 

2280 

2290 

230C 

2310 

2320 

2330 

2370 

2380 

2390 

2400 

2410 

2420 

2430 

2440 

2450 

2460 

2470 

2480 

2490 

2500 

2510 

2520 

2530 

2540 

2550 

2560 

2570 

2530 

2590 

2520 

2510 

2620 

2630 

2640 

2650 

2660 

2670 

2680 

2690 

2700 

2710 


FOR  I>1  TO  Ibtuif 
Kttp(2 , I)*Icol 
R«<ly«R««ly'»'Add(  1 ) 

Add(t)-Addtl) 

IF  R*«ly>Irow+l  THEN  Irow»Irow+l 
Addy ( I ) »R«« ly- I row 
K*ep(l,I)-Irow 
Icol»Icol+l 
NEXT  I 


*****SET  UP  INDEX  OF  REFRACTION  FOR  PATH  OF  INCIDENT  BEAM  **** 


SET  UP  CRYSTAL  TO  MEAN  INDEX 

M#an»ABS ( ( I top+Ibot ) /2 ) 

FOR  I«1  TO  I size 

FOR  J«1  TO  I size 

Icrys(I , J)»Mean 
Itrac(I, J)*0 
NEXT  J 
NEXT  I 
! 

! 

! 

Grad«(Itop-Ibot)/2  I  INDEX  OF  REFRACTION  GRADIENT 
FOR  1*1  TO  Ibeamf 
Irow=Keep(l,I) 

Icol*Keep(2,I) 

Crys*Grad*SIN(((Ibeamf/2)-(I-l))*(PI/Ibeamf)) 
Crys*Crys+Mean 
Icrys(Irow,Icol)»Crys 
NEXT  I 
I 
! 
j 

i  SET  UP  VALUES  DOWN  THE  BEAM  PATH 
! 

FOR  1=1  TO  Ibeamf 
Icol*Keep(2,I) 

Irow=Keep(l,I) 

Imdex=Icrys(Irow,Icol) 

Imdexp(I)=Imdex 

Itemp=Icol 

Realx=0 

FOR  J»Irow  TO  Isize 

I c  ry s ( J , I temp ) =Imdex 
Itrac(J,Itemp)=l 
Realx=RealX't'Add(  I ) 

Itemp=Icol-Realx 


•lEx: 


'.rEi-rr 


IF  Itemp>Isi2e  THEN  Itemp 
IF  I’:emD<0  THEN  ItemD=0 


=Isi2e 


SET  THE  STARTING  COLUMNS  FOR  THE  WAVEFRONT  AND 
INITIALIZE  FLAGS 


FOR  1=1  TO  IBEAMF 

Ipermc(I)=Keep(2,I) 

lmarkba)*0 


2720 

2730 

2740 

2750 

2760 

2770 

2780 

2790 

2800 

2810 

2820 

2830 

2840 

2850 

2860 

2870 

2880 

2890 

2900 

2910 

2920 

2930 

2940 

2950 

2960 

2970 

2980 

2990 

3000 

3010 

3020 

3030 

3040 

3050 

3060 

3070 

3080 

3090 

3100 

3110 

3120 

3130 

3140 

3150 

3160 

3170 

3180 

3190 

3200 

3210 

3220 

3230 

3240 

3250 

3260 

3270 

3280 

3290 

3300 

3310 

3320 

3330 

3340 

3350 

3360 

3370 

3380 

3390 

3400 

3410 


lm«rki(l)«0 
NIXT  I 


HOVEKENT  ^******^aaaaaaaaaaaaaaaaaaaaaaaaaaaa 

K4XE  SURE  BEAM  IS  STILL  IN  CRYSTAL 

FOR  !■!  TO  Ibcamf 

IF  Im«rkb(I)-l  THEN 

IF  Keep (1,1) <2  THEN  GOTO  5140 
END  IF 

IF  Imerks(I)«l  THEN 

IF  Keep(2,I}>Islze-10  THEN  GOTO  5140 
END  IF 
NEXT  I 
1 

I  aa^’^aa^HOVEMENT  aaa^a***a*a*aaaaaaaaaaaaaaaaaaaaaaaaaaaaaa 

1 

t 

WHILE  Rcheck<Rstop 
1 

FOR  I>1  TO  Ibeamf 
Icol«Keep(2,I) 

Irow«Keep(l jl) 

Imdex»Icrys ( I row , Icol) 

! 

!  HARK  REFLECTION  IF  NECESSARY 
1 

IF  Icol<=l  THEN 

IF  Imarks(I)=0  THEN 
Imarks(I)»l 
Realad(I)»l 
END  IF 
END  IF 

IF  Irow  =Isi2e  THEN 
Imarkb(I)=l 
END  IF 

! 

! 

I  MOVE  IF  IT  IS  TIME 
! 

IF  Itime(I)<Imdex  THEN 
ItimevI)=Itirae(I)+l 
ELSE 

Itime(I)=0 

Reaiaa( I )=Realad(I ) +Add( I ) 

IF  Imarks(I)=0  TlteN 

Icol=Ipennc(I)-Realad(I) 

IF  lcol<0  THEN  lcol=0 
ELSE 

Icoi=Reaiad(I) 

END  :f 

Keep^2,I0=Icol 
IF  Imarkab(I)=0  THEN 
Keep(l,I)=Irow+l 

Sep(l,I)=Irow-l 
END  IF 


END  IF 
NEXT  I 


3420  I 

3430  I  PUT  POINTS  IN  FILE  AND  DISPLAY  GRAPHICS 

3440  I  IF  IT  IS  TIME  TO  DO  SO 

34S0  I 
3460  I 

3470  Igraph>Igraph<«'l 
3480  IF  Igrapn>Gnnax  THEN 
3490  Icout«lcout+l 

3500  FOR  L»1  TO  Ibeamf 

3510  lgraph«0 

3520  Iout(0,L,Icout)«Keep(2.L)  !X  COORD  OF  POINT 

3530  Iout(l,L,Icout)=Isize-Keep(l.L)  !Y  COORD  OF  POINT 

3540  PENUP 

3550  PLOT  Iout(0,L,Icout) ,Iout(l,L,Icout) 

3560  NEXT  L 

3570  END  IF 
3580  I 
3590  I 
3600  I 
3610  I 
3620  I 

3630  Rcheck^Rcheck-*-! 

3640  t 

3650  !  CHECK  IF  TIME  TO  RECOMPUTE  DIRECTION  ****************** 

3660  I 

3670  NDIR*Mdir+l 
3680  IF  Mdir>«Idir  THEN 
3690  Mdir=0 

3700  PRINT  TABXY(1,11),  CURRENT  TIME  PERIOD  »  ,  Rcheck 

3710  ! 

3720  I 
3730  ! 

3740  1 

3750  1 

3760  1 

3770  I 

3780  FOR  K*0  TO  Piece-1 
3790  ! 

3800  i  ****************  SET  FOR  5  PIECES  ***************** 

3810  I 

3820  IF  K»0  THEN  KkaO 

3830  IF  K*1  THEN  Kk=4 

3840  IF  R>2  THEN  Kk^l 

3350  IF  Ka3  THEN  Kk=3 

3860  IF  K«4  THEN  Kka2 

3870  !IF  K=5  THEN  Kk=7 

3830  I  IF  K=6  THEN  Kka3 

3890  !IF  Ka7  THEN  Kk=6 

3900  I  IF  K=8  THEN  Kka4 

3910  !IF  K=9  THEN  Kk«5 

3920  This=Kk*Opl+l 

3930  Thata(Kk+l)*Opl 

3940  IF  Kk»Piece-l  THEN  That«Ibeamf 

3950  I 

3960  ! 

3970  J  FIRST  CHECK  TO  MAKE  SURE  THE  FRONT  IS  NOT  PARTIALLY  REFLECTED 
3930  ! 

3990  ; 

4000  lsxipb=0 

4010  I  jev;.-nark=0 

4020  IF  Reflchkb(K)=0  THEN 

4030  FOR  M=This  TO  That 

4040  IF  Imarkb(M)*l  THEN 

4050  Isetmarkal 

4060  FOR  N=This  TO  That 

4070  IF  Imarkb(N)*0  THEN  Iskipb*! 

4080  NEXT  N 

4090  END  IF 

4100  NEXT  M 

4110  END  IF 


4120 
4130 
4140 
4150 
4160 
4170 
4180 
4190 
4200 
4210 
4220 
4230 
4240 
4250 
4260 
4270 
4230 
4290 
4300 
4310 
4320 
4330 
4340 
4350 
4360 
4370 
4380 
4390 
4400 
4410 
4420 
4430 
4440 
4450 
4460 
4470 
4480 
4490 
4500 
4510 
4520 
4530 
4540 
4550 
4560 
4570 
4530 
4590 
4600 
4610 
4620 
4630 
4640 
4650 
4660 
4670 
4630 
'  ' 
T  J  - 

'k  I 

4720 

4730 

4740 

4750 

4760 

4770 

4780 

4790 

4800 

4810 


IF  Isetnark>l  AND  lskipb«0  THEN  Rtflchkb(K)«l 

lskips»0 

lsetin«rk>0 

IF  Reflchki(K)-0  THEN 
FOR  M»This  TO  That 
IF  Imarks7H)«l  THEN 
FOR  N»Thia  TO  That 
Isatmark^l 

IF  l!narks(N)»0  THEN  Iskips*! 

NEXT  N 
END  IF 
NEXT  M 


IF  Isetmark-l  AND  lsklps«0  THEN  Reflchks(K)«l!SKIP  THIS  IF  NOT  NECESSARY 


END  IF 


IF  lakips«0  AND  lakipb«0  THEN  IIF  WAVEFRONT  IS  NOT  PARTIALLY  REFLECTri 


SLOPE  OF  PIECE  OF  WAVEFRONT  IS  YY/XX 


Yy*^Ke«p(l ,This)+Addy(This) )“(Keep(l , That|+Addy(That) ) 

71* .  . -  .  .  - 


I 


|t^e{Tjiis)/^Icrys|kaep^i  i 
THEN 


Y2*I tima (That )/Icrys (Keep 
Yj*Yy+Yl-Y2 


This ) , Keep ( 2 , This ) ) 
That) , Keep (2, That) ) 


Imarks(This)aO 

Xx» ( Ipermc (That ) -Realad(That ) )  -  ( Ipertnc (This )  -Realad(This ) ) 

ELSE 

Xx*'Realad(That )  >Realad(This ) 

END  IF 

IF  Xx>l  THEN 
FOR  G»This  TO  That 

Add(G)>ABS(Yy/Xx)  I  FIGURE  NEW  DIRECTION 
Caxis(G)«Xx/Opl 

IF  Kk»Piece-l  THEN  Caxis(G)»Xx/ (Ibeamf-(K*Opl ) ) 

NEXT  G 

PRINT  Add(G-l) 

END  IF 


t 


! 


FOR  I«This  TO  That 
Icol»Keep(2,I) 
Irow»Keep(l , I) 
Imdex»Imdexp(I ) 


i\ea-.k=J 


Oi. 


I  CHECK  FOR  REFLECTION 


I 


IF  Imarkb(I)»0  THEN 

FOR  J*Irow  TO  Isize 

Icrys( J,Itemp)«Imdex 
Realx*RealX'«'Add(  I ) 

I tempalcol'Realx 
IF  ltemp<0  THEN  ltemp«0 


3 


IF  Ittmp>Isiz«  THEN  It«inp«Isize 
Itr«c(J,Itemp)*l 
NEXT  J 


! INCIDENT  BEAN  DOMINATES 
FOR  J»0  TO  I  row 

Ichng*Irow-J 

IF  Itrac( Ichang, ltemp)»0  THEN 
Icrys{lchng, Itenm;»Imdex 
ELSE 

Icrys(Ichng,  Iten>p)»Icrys(Ichng,  I  --itip) 
END  IF 
t 

Realx«Realx-Add( I ) 

Itempsicol-Raalx 
IF  ltemp<0  THEN  lteinp«0 
IF  Itemp>Isize  THEN  itemp*Isize 
NEXT  J 


END  IF 
NEXT  I 
END  IF 
NEXT  K 
END  IF 


I  CHECK  FOR  END  OF  RUN  >******>************** 

END  WHILE 

T2-TIMEDATE 

PRINT  TOTAL  RUN  TIME  *  ,T2-T1 
GCSUB  5360  1  WRITE  FILE  TO  DISK 

STOP 

!  SUBROUTINE  TO  INITIALIZE  GRAPHICS 

GRAPHICS  ON 

GCLEAR 

SHOW  0,lsiz«,0,lsiza 
PLOT  0,0, -2 
FOR  Q-6  to  I size 
PLOT  0,0, -I 
NEXT  Q 

FOR  Z»0  TO  I size 
PLOT  Isize,Z,-l 
NEXT  Z 

FOR  Q»Isize  TO  0  STEP  -I 
PLOT  Q,  Isize,-1 
NEXT  Q 

FOR  Z«Isize  TO  0  STEP  -1 

PLOT  0,2, -1 

NEXT  2 

PLOT  0,0, -2 

RETURN 

!  SUBROUTINE  TO  OUTPUT  FILE  TO  DISK 
MASS  STORAGE  IS  :,  700,0 
CREATE  BOAT  NameS  10404,2 
A55I0N  Pa:.'.  1  TO  'larrsO 


!  SUBROUTINE  TO  OUTPUT  FILE  IF  ERROR  OCCURS 

MASS  STORAGE  IS  700,0 

CREATE  BOAT  NameS , 10404 , 2 

ASSIGN  Path  1  TO  NameS 

OUTPUT  Path  l;Iout(’») 

ASSIGN  Path  1  TO  * 

PRINT  error!  PROGRAM  TERMINATED 
T2-TINEDATE 

PRINT  TOTAL  RUN  TIME  •  .T2-T1 
STOP 


'  1/ 
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